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Abstract 



X 

$—i ' The possibility of additional quarks and leptons beyond the three gener- 



ations already established is discussed. The make-up of this Report is (I) 
Introduction: the motivations for believing that the present litany of elemen- 
tary fermions is not complete; (II) Quantum Numbers: possible assignments 
for additional fermions; (III) Masses and Mixing Angles: mass limits from 
precision electroweak data, vacuum stability and perturbative gauge unifica- 
tion; empirical constraints on mixing angles; (IV) Lifetimes and Decay Modes: 
their dependence on the mass spectrum and mixing angles of the additional 
quarks and leptons; the possibility of exceptionally long lifetimes; (V) Dy- 
namical Symmetry Breaking: the significance of the top quark and other 
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heavy fermions for alternatives to the elementary Higgs Boson; (VI) CP Vi- 
olation: extensions to more generations and how strong CP may be solved 
by additional quarks; (VII) Experimental Searches: present status and future 
prospects; (VIII) Conclusions. 
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I. INTRODUCTION. 



The elementary spin-half fermions as we now know them are the quarks and leptons. 
The principal constituents of normal atoms and normal matter are the electron, as well as 
the up and down quarks which comprise the valence quarks of the protons and neutrons in 
the atomic nucleus. In addition, there is the electron neutrino which was first postulated 
by Pauli M in 1931 to explain conservation of energy and angular momentum in nuclear 
/3-decay, and which was eventually discovered in 1956 by Reines and Cowan 0. 

These quarks and leptons - the up and down quarks, the electron and its neutrino - 
comprise what is now called the first generation. The first intimation that Nature is more 
complex came with the discovery of the muon in 1937 P,f§]. The muon appears identical 
to the electron except for its mass which is ~ 200 times heavier. The muon appeared so 
surprising that there was a famous comment by I.I. Rabi ||: "Who ordered that?" 

The fact that the muon neutrino differs from the electron neutrino was established in 
1962 ||. The strange quark had already been discovered implicitly through the discovery of 
strange particles beginning in 1944 PH9|. Completion of the second family with the charmed 



quark, predicted in 1970 [10|, was accomplished in 1974 [p|-|T4|| . At first only hidden charm 
was accessible but two years later explicit charm was detected ||15|| . 

By this time, a renormalizable gauge theory was available |T6T-|2"0"[ based on the first two 



generations and incorporating the Cabibbo mixing J2TJ] between the two generations. 

The situation became even more challenging to theorists when experimental discovery 
of the third generation of quarks and leptons began with the tau lepton, discovered in 1975 



22[ in e + e scattering at SLAC. Next was the bottom quark in 1977 p3| , |24| . The top quark, 



at ~ 175 GeV much heavier than originally expected, was finally discovered in 1995 |25| , |26l . 
Together with the r neutrino which presumably participates (its distinct identity - while not 
questioned - is not fully demonstrated) in tau decay, this completed the third generation. 

Since the present review is dedicated to the premise of further quarks and leptons beyond 
the third generation, it is worthwhile to recall to what extent and how the third generation 
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was anticipated from the existence of the first two generations, why it is regarded as the end 
of the litany of quarks and leptons and what loopholes there are in the latter arguments. 

One early theoretical anticipation of a third generation was the paper of Kobayashi and 
Maskawa [^7] who pointed out, at a time (1973) when only three flavors u, d, s of quark were 
established, that the existence of six flavors in three generations would allow the standard 
model naturally to accommodate CP violation. 

Study of the formation of the light elements (Hydrogen, Deuterium, Helium, and 
Lithium) in the early universe was started earlier in the 1960's f||[29l. In the 1970's tighter 
constraints were found based on the steadily-improving estimates of the primordial abun- 
dances of these light isotopes. Since the expansion rate of the universe in this era of Big-Bang 
Nucleosynthesis, and hence the abundances, depends sensitively on the number of light neu- 
trinos it was then possible to limit the acceptable number. The group of Schramm et al. 
|30H32| found in this way that the number of generations should not be greater than four 
pOfl , or in some analyses not greater than three (see e.g. footnote 4 on page 242 of |]3T|); it is 
surely remarkable that such a strong constraint was found from early universe considerations 
already in 1979, a decade before the situation was clarified using colliders. 

A current plot of the primordial A He abundance (whose exact value is still controversial 
in 1998) versus neutrino number n(u) for n{v) = 3.0, 3.2, 3.4 is given in Figure 1. The main 
point is that the neutrino number from cosmology is by now tied very closely to the high- 
energy experimental value in what is the strongest known link between particle theory and 
cosmology. 
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FIG. 1. Helium-4 production for N v = 3.0,3.2,3.4. The vertical band indicates the baryon 
density consistent with (D/H)p = (2.7 ± 0.6) x 10 -5 and the horizontal line indicates a primeval 
Helium-4 abundance of 25%. The widths of the curves indicate the two-sigma theoretical uncer- 



tainty. Figure from Ref. |33 ]. 



In 1989, there came an experimental epiphany concerning the number of generations, 
or more precisely, the number of light neutrinos. This arose from the measurement of the 



Z width at SLAC |34| , |35|| and especially at CERN p6|-|39||. The answer from this source is 
indisputably equal to three. The argument is simple: One can measure the total width 
of the Z to high accuracy, and then subtract the visible width to get the invisible width. 
Identifying the invisible width with neutrino decays leads to [STT: 
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n(v) = 3.00 ± 0.025 (1) 

This provides compelling proof that there are only three conventional neutrinos with mass 
below Mz/2 ~ 45 GeV. And, by extrapolation, it leads to the idea that there are only three 
quark-lepton families. Since this ties in nicely with the KM mechanism of CP violation and 
with the Big-Bang-Nucleosynthesis, the overall picture looks very attractive. 



However, this finality was not universally accepted [0,0- There are other reasons for 



entertaining the possibility of further quarks and leptons: 

• In many grand-unified theories (GUTs), there naturally occur additional fermions. 
Although the minimal SU(5) GUT can contain only the basic three families, extension 
to any higher GUT such as SO (10) or E(6) inevitably adds new fermions. In 5*0(10) 
this may be only right-handed neutrinos (leptons) but in E(6) there are also non-chiral 
color triplets(quarks) and color singlet (leptons). Although there is no direct evidence 
for GUTs they are attractive theoretically and suggestive of how the standard nodel 
is extended. 

• Models of CP violation which solve the strong CP problem without axions generically 
require additional quarks. 



For some models of dynamical symmetry breaking [[43H46| the top quark mass, al- 
though higher than originally expected, is still not quite high enough to play its role 
in electroweak symmetry breaking. This might be interpreted as evidence for even 
heavier quarks. 

It has been shown that a non-supersymmetric model with four generations can have 
successful unification of the gauge couplings at the unification scale. 

In recently-popular models of gauge-mediated supersymmetry breaking, additional 
vectorlike quarks and leptons arise automatically. In addition, in models in which 
higher dimensions arise at the TeV scale [47], it has been shown that if some 
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standard model fields live in the higher dimensional space, low-scale gauge unification 
can be obtained-the Kaluza-Klein excitations of these fields, which could be rather 
light, must be vectorlike. 

None of these reasons is fully compelling but each is suggestive that one should keep 
alive the study of this issue. Our hope is that this review will play a role in encouraging 
further thought about this open question. 

The present review contains the following subsections: Section II is on the possible quan- 
tum numbers of additional quarks and leptons; Section III discusses their masses and mixing 
angles; Section IV deals with lifetimes and decay modes; Dynamical symmetry breaking is 
in Section V; CP violation is treated in Section VI, and finally in Section VII there is a 
treatment of the experimental situation and in Section VIII are the conclusions. 
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II. QUANTUM NUMBERS. 



When we add fermions to the standard model, there are choices in the possible quantum 
numbers. 

Under the color 577(3) c group we refer to color triplets as quarks, color antitriplets as 
antiquarks. Color singlets which do not experience the strong interaction are generically 
referred to as leptons. Higher representations of color such as 6, 6, 8,... may be called 
quixes, antiquixes, queights, and so on. Such exotic color states are necessary in some 



models to cancel chiral anomalies. For example, in chiral color |49|j50|] one version (called 
Mark II in |49[]) involves three conventional fermion generations, an extra Q = 2/3 quark, 
and an SU(3) sextet fermion or quix. The extended gauge group of chiral color is 577(3)/, X 
SU(3) R x SU(2) L x U{1) Y and we may list the fermions by their (SU(3) L , SU(3) R , Q) 
quantum numbers. There are three colored weak doublets 



3[(3,1,|) + (3,1,-1)], (2 ) 



eight colored weak doublets 



a weak singlet quix 



4(1, 3,-|) + 3(1,3,+!) + (3,1,|), (3) 



(6, l,-^) + (l, 6, i), (4) 



and three charged leptons and their neutrinos. The quix plays an essential role in anomaly 
cancellation. But in this review we restrict our attention only to quarks and leptons because 
while more exotic color states are a logical possibility it is one which is difficult to categorize 
systematically. 

Quarks and leptons may be either chiral or non-chiral. The latter are sometimes alter- 
natively called vector-like. Let us therefore define the meaning of these adjectives. 

Chiral fermions are, for present purposes, spin-| paricles where the left and right com- 
ponents transform differently under the electroweak gauge group SU(2)l x U(l). All the 
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fermions of the standard model are chiral. This means that they are strictly massless before 
the electroweak symmetry is broken. 

The simplest generalization of the standard model is surely to add a fourth sequential 
family. Of course, a fourth light neutrino is an immediate phenomenological problem with 
the invisible Z width, but the addition of a right-handed neutrino can resolve this. 

More generally we may add a chiral doublet quark or lepton where the left-handed 
components transform as a doublet of SU(2) L and the right-handed components as singlets. 
A chiral doublet of quarks is: 



U R ; D R 



while a chiral doublet of leptons is 



N, 



R, 



Equally possible are chiral singlets such as 



E 



R 



(5) 



(6) 



U T 



(7) 



or 



D, 



(8) 



for quarks, or for leptons 



Nr 



(9) 



or 



(10) 



Of course, with chiral doublets or singlets, there is a constant danger of chiral anomalies. 
In the standard model, there is a spectacular cancellation in each generation between the 
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anomalies of the chiral doublets of quarks and leptons. In other models one sometimes adds 
mirror chiral doublets to cancel anomalies. For example, a mirror chiral doublet of quarks is 



(11) 



There can also be non-chiral (also known as vector-like) fermions, where the right and 
left components transform similarly under the electroweak SU(2) x U (1) group. For example 
a vector-like quark doublet is 



(12) 











1 

L 





R 



A vector-like doublet is not to be confused with the doublets occurring in the left-right 
model |5TH53| where the gauge group is extended to SU(2)l X SU(2)r x U{1)b-l- The 
U and D quarks transform there chirally as in (|^) above under SU(2)l, but the right- 
handed singlets Ur and Dr of SU{2)l are then assumed to transform as a doublet under 
the additional gauge group SU{2) R . 

General patterns for adding anomaly-free charge-vectorial chiral sets of fermions which 
acquire mass by coupling to the Higgs doublet of the standard model have been studied in 



54J and further developed in . An extensive analysis of the possible quantum numbers 

They looked at the general structure of 



for additional fermions can be found in Ref. 
exotic generations given the gauge and Higgs structure of the standard model. A similar 



analysis for left-right symmetric models was carried out in Ref. |58 . 

In grand unified theories (GUTs) all types of additional fermions are possible. For exam- 
ple, in SU(5) there may be no additional fermions. But in 5*0(10) there must be at least an 
additional chiral right-handed neutrino in each family. In E(6) each 27 of fermions contains 
not only two extra neutrino-like states but also a 5 + 5 of SU (5) which contains a non-chiral 
singet of quarks and a non-chiral doublet of leptons. 

In superstring models, and M-theory models, the additional fermions are even less con- 
strained. For example E(G) and its content has been a familiar superstring possibility since 
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the beginning ||, and any gauge group contained in at least 0(44) is possible. Families and 
anti (mirror) families often occur in superstrings. 

Some extensions of the standard model require extra chiral fermions to cancel anomalies 

e.g. 



As already mentioned, chiral color [|5], anomaly cancellation dictates addition of 
quixes. 



• In some GUTs, e.g. SU(15) anomaly cancellation requires |Kj inclusion of mirror 
fermions. 

More interesting are the types of additional quarks and leptons appearing in extensions of 
the standard model motivated by attempting to explain shortcomings of the model itself. 
Examples are: 



Attempts to solve the strong CP problem without an axion pTj , ^ can introduce a 
non-chiral doublet of quarks, or non-chiral singlets (but not both), 

In trying to explain the three families through anomaly cancellation, the 331-model 
[]63| , |64ll extends the individual families of the standard model by adding non-chiral 
singlets of quarks. The charges of the additional quarks differ between the families. 
In a sense, this is not "Quarks and Leptons beyond the Third Generation" since the 
quarks are being added to the discovered generations. Nevertheless, our framework is 
sufficiently general to accommodate this possibility - as additional non-chiral singlets 
of quarks. 

In the non-supersymmetric standard model, gauge unification of the couplings fail-the 
three couplings do not meet at a point. It has been noted |6l| that extending the model 
to allow a fourth generation introduces enough flexibility that a successful unification 
of couplings can occur. 
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• A potential problem of the minimal supersymmetric standard model is that the scalar 
quark masses must be very nearly degenerate to avoid large tree-level flavor-changing 
neutral currents. This degeneracy is very unnatural, using the conventional mecha- 
nism of gravitationally mediated supersymmetry breaking (where the supersymmetry 
breaking part of the Lagrangian is transmitted to the known sector via gravitational in- 
teractions). However, in gauge- mediated supersymmetry breaking, the supersymmetry 
breaking is transmitted to the known sector via the gauge interactions of "messenger" 
fields. Since the gauge interactions are flavor-blind, the degeneracy of squark masses 
is natural. The messenger fields are, in the simplest case, composed of a 5 + 5 of 
SU(5) (or of several such fields), which constitute a vectorlike Q = —1/3 quark and a 
vectorlike lepton. 

• There has been recent excitement about the possibility that additional spacetime di- 
mensions could be compactified at (or even well below) the TeV scale ||47|| . In such 
models, the Kaluza-Klein excitations must be vectorlike (to avoid having, for example, 
too many light neutrinos). 

In this Report, we will primarily concentrate on chiral quarks and leptons, or vector like 
quarks and leptons, since they appear in the majority of models. However, when appropriate, 
we will note how our various constraints and bounds will apply to mirror quarks. 
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III. MASSES AND MIXING ANGLES. 



One of the most unsatisfactory features of the standard model is the apparent arbitrari- 
ness of the masses and mixing angles of the known fermions. Although masses and mixing 
angles can be accommodated in the standard model (with addition of right handed neutrino 
fields, if necessary), there is no understanding of their values. An entire industry of model- 
building has developed in an attempt to provide some theoretical guidance, ranging from 
flavor symmetries to relationships from grand unification, but no model seems particularly 
compelling. 

In the case of additional fermions, the masses and mixing angles also are arbitrary. 
Nonetheless, some general features can be found. Phenomenological bounds can be obtained 
from high precision electroweak studies, theoretical bounds can be obtained from requiring 
the stability of the standard model vacuum and from requiring that perturbation theory be 
valid (up to some scale). In this Section, these bounds are explored in some detail. We 
will start with a discussion of the phenomenological constraints from precision electroweak 
studies. Then we will consider bounds from vacuum stability and perturbation theory. In 
the vast majority of analyses of these bounds, the authors focused on bounds to the top 
quark mass, since its mass was unknown until relatively recently, thus we will first look at 
constraints on the top quark mass, and then generalize the results to find bounds on masses 
of additional quarks and leptons. Finally, we will discuss plausible models for the mixing 
angles. 



A. Precision Electroweak Constraints 



1. Chiral Fermions 



In the past decade, high precision electroweak measurements have led to remarkable 
constraints on potential physics beyond the standard model. The most important of these is 
the p parameter. As originally pointed out by Veltman p6|,|67| , the tree level mass relation 
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is very sensitive to non-standard model physics (ruling out, for example, significant vacuum 
expectation values for Higgs triplets). Since the relation is good to better than 1%, it is 
assumed that deviations from the relation are due to electroweak radiative corrections, which 
are sensitive to new particles in loops. 

An extensive, detailed analysis of electroweak radiative corrections can be found (with a 
long list of references) in the work of Peskin and Takeuchi |}8],|69| . They define the S and T 
parameters as 

aS =4e 2 [n'33(0)-n'3 Q (0)] 

« T -^£M-[ n n(o)-n 33 (o)] (14) 

where a is the fine-structure constant, Hxy{q 2 ) is the vacuum-polarization amplitude with 
(XY) = (11), (22), (33), (3Q), (QQ), and IT = n(g2) g ~ n(0) . Roughly, T is a measure of the 
deviation of the p parameter from unity, coming from isospin violating contributions. S is 
an isospin symmetric quantity which measures the momentum dependence of the vacuum 
polarization; it is roughly the "size" of the new physics sector[] 

As an example, Peskin and Takeuchi consider the case of a chiral lepton doublet, E and 
N. In the limit Mn,Me » Mz, they show that if the mass splitting is small, then the 
contribution for S is just tt-, and the contribution for T is — \ A si K,i . For a chiral quark 

J 07r' 127rsuo QwMyy ^ 

doublet, these contributions are tripled. As stated above, one can see that T is a measure of 
the isospin splitting, while S is a measure of the size of the new sector. Thus, for a complete 
degenerate fourth generation, the contribution to S is 

S = ^--0.21 (15) 
on 

while the contribution to T is 



x An alternative representation, using parameters ex, £2 and 63 can be found in Ref. [70]. 
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-I ! «J (IQ) 

12tt sin 2 6 w M$v 47rsin 2 # Ty M^ v ; 

It is important to note that these results were obtained under the assumption that the extra 
fermion masses were much greater than that of the Z, and that the mass splitting was small. 
Peskin and Takeuchi give more exact expressions. 

What are the experimental bounds? The values can be determined from the data on 
Z-pole asymmetries, M z , M w , R h m top , R b and Z-decay widths. The contribution from 
new physics, S new = S — Ssm an d T new = T — Tsm, can be determined [|4*0] , |n] , [72]l . The 



value of this contribution depends on the Higgs and top quark masses (which affect Ssm, 
for example). For a top mass of 175 GeV and a Higgs mass between Mz and 1 TeV, Erler 
and Langacker W% find 



S new = -0.20i°;*° (3a) (17) 

One can immediately see a apparent conflict between Eqs. [15] and [17|. Independent of the 
mass, a fourth generation chiral multiplet is roughly 3 standard deviations off. This led Erler 
and Langacker to claim that a fourth sequential family is excluded at the 99.8% confidence 
level. 

However, there are several reasons why we believe it is premature to exclude a fourth 
sequential family. First, of course, is the fact that many 3a effects in recent years have 
disappeared. Second, the results for S are based on virtual heavy fermion loops. Any 
additional new physics will likely make a similar contribution. For example, Erler and 
Langacker show that the allowed range for S in minimal SUSY is —0.171^2, where the 
error bars are la, and this is only in conflict by 2.2 standard deviations. Third, it has been 
noted |73| that Majorana neutrino masses (which may be needed to give neutrino masses in 



the right range) lower S, as do models which involve mixing of two scalar doublets [[F4| , and 
this could reduce the discrepancy further. Finally, the result above assumed a degenerate 
family. If one uses the exact expressions, takes M H = 180 GeV, M D = M v = 150 GeV, 
M/v = 100 GeV and Me = 200 GeV, for example, one finds the contribution to S to be 
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approximately 0.11, rather than 0.21. This also would lower the discrepancy to slightly more 
than 2a. 

For S ~ 0.2, the 2a upper bound on T is approximately 0.2, leading to bounds on the 
mass splitting from Eq. [TJ. For quarks, the splitting must be less than 54(Mw/Md) 
GeV and for leptons must be less than 162 {My// Me) GeV, at the 2a level. Note that for 
quark masses above 200 GeV, the splitting must be less than 10 percent. 

Of course, the extra fermions must not contribute significantly to the width of the Z, 
and their masses are thus bounded from below by Mz/2. Other experimental bounds will 
be discussed in Section VII. 

2. Non-chiral Doublets 

Vector-like fermions do not contribute in leading order to 5* and T, and thus the values 
of these parameters do not constrain their masses. However, since vector-like doublets do 
not couple to the Higgs boson, the mass terms involving the E and the N cannot violate 
isospin invariance, and thus the masses must be degenerate at tree-level, as must the masses 
of the U and D. Even if one adds a singlet Higgs field, the degeneracy will remain. Only 
a Higgs triplet can split this degeneracy, but Higgs triplet vacuum expectation values are 
severely constrained by the p parameter. We conclude (in the absence of sizeable mixing 
with lighter generations) that the masses of states in a vector-like doublet are degenerate 
at tree level. The masses will be split by a few hundred MeV due to electroweak radiative 
corrections-this calculation will be done in the next Section. 

3. Other Fermions 

Non-chiral singlets will have arbitrary mass terms and arbitrary couplings to any Higgs 
singlets. No constraints can be placed on their masses. 
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B. Vacuum Stability Bounds 



Upper bounds on fermion masses can be obtained from the requirement that fermionic 
corrections to the effective potential not destabilize the standard model vacuum. We will 
first discuss the effective potential, and its renormalization group improvement. Then, the 
bounds from the requirement of vacuum stability will be discussed, first for the top quark 
mass, and then for additional quarks and leptons. 

1. The Effective Potential 

An extensive review of the effective potential and bounds from vacuum stability ap- 
peared in 1989 [[75]. Since then, the potential has been improved, including a proper 
renormalization-group improvement of scalar loops, and the bounds have been refined to 
much higher precision. In addition, the discovery of the top quark has narrowed the re- 
gion of parameter space that must be considered. In this section, we discuss the effective 
potential and its renormalization-group improvement. 

It is easy to see how bounds on fermion masses can arise. The one-loop effective potential, 
as originally written down by Coleman and Weinberg |76|] can be written, in the direction 
of the physical field, as 

V{4>) = V + V l (18) 

where 

Vo = + ~A0* (19) 

and 

*=BS5E<-l>W(*Jta^ < 20 > 

i 

where the sum is over all particles in the model, F is the fermion number, rji is the number of 
degrees of freedom of the field i, and Ai 2 ((p c ) is the mass that the field has in the vacuum in 

17 



which the scalar field has a value C . In the expression for Vi, we have ignored terms which 
can be absorbed into Vo-these will be fixed by the renormalization procedure. In the standard 
model, we have for the W-boson, M 2 (<p c ) = \g 2 <Pl, for the Z-boson, M 2 (4> c ) = \(g 2 + g /2 )0 2 , 
for the Higgs boson, M 2 {4> c ) = — /^ 2 + 3A0 2 , for the Goldstone bosons, M 2 (<p c ) = — /^ 2 + A0 2 
and for the top quark Ai 2 (4> c ) = \h 2 <t> 2 - For a very large values of 0, quadratic terms are 
negligible and the potential becomes 

V = ^A0 4 + B(f) 2 ln(0 2 /M 2 ) (21) 

where 

B = ^[4A 2 + 1(3^ 4 + 2g 2 g' 2 + g") - h*} (22) 

One can see that if the top quark is very heavy, then h is large and thus B is negative. In 
this case, the potential is unbounded from below at large values of 0. This is the origin of 
the instability of the vacuum caused by a heavy quark. 

Although this form of the effective potential is well known, it is NOT useful in determin- 
ing vacuum stability bounds. The reason is as follows. Suppose one denotes the largest of 
the couplings in a theory by a, in the standard model,for example, a = [max(A, g 2 , h 2 )]/ (47r). 
The loop expansion is an expansion in powers of a, but is also an expansion in powers of 
logarithms of (f) 2 /M 2 , since each momentum integration can contain a single logarithmic 
divergence, which turns into a ln(0 2 /M 2 ) upon renormalization. Thus the n-loop potential 
will have terms of order 

a n+1 [ln(0 2 /M 2 )] n (23) 

In order for the loop expansion to be reliable, the expansion parameter must be smaller 
than one. M can be chosen to make the logarithm small for any specific value of the field, 
but if one is interested in the potential over a range from <pi to 02, then it is necessary for 
aln(0 1 /0 2 ) to be smaller than one. In examining vacuum stability, one must look at the 
potential at very large scales, as well as the electroweak scale, and the logarithm is generally 
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quite large. Thus, any results obtained from the loop expansion are unreliable; in fact, the 
bound on the top quark mass can be off by more than a factor of two. 

A better expansion, which does not have large logarithms, comes from solving the renor- 
malization group equation (RGE) for the effective potential. This equation is nothing other 
than the statement that the potential cannot be affected by a change in the arbitrary pa- 
rameter, M, i.e. dV/dM = 0. Using the chain rule, this is 

wm + ^wr^ v = < 24) 

where (5 = Mdgi/dM and there is a beta function for every coupling and mass term in the 
theory. The 7 function is the anomalous dimension. 

It is important to note that the renormalization group equation is exact and no approxi- 
mations have been made. If one knew the beta functions and anomalous dimensions exactly, 
one could solve the RGE exactly and determine the full potential at all scales. Although we 
do not know the exact beta functions and anomalous dimensions, we do have expressions 
for them as expansions in couplings. Thus, by only assuming that the couplings are small, 
the beta functions and 7 can be determined to any level of accuracy and V((f>) can be found. 
The resulting potential will be accurate if gi « 1 and will not require ^ln(0/M) << I. 

For example, in massless A0 4 theory, the RGE can be solved exactly to give 

V = iA'(t,A)G 4 (t,A)0 4 (25) 
where t = ln(0/M) and X'(t, A) is defined to be the solution of the equation 

dt - (l + 7(A')) {b) 
with the boundary condition being determined by the renormalization condition. G(t, A) is 
defined as exp(— 4 J* * dt'('y(X')/(l + 7(A'))). Note that this potential gives the same result as 
before in the limit that 7 = and (3 = constant. Then G — 1 and A' = f3t + constant. With 
t = ln(0/M) this gives the 4 ln(0/M) terms as above. 
What about the massive case? The RGE is given by 
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^ M m + "4 + + >^W' ~ 7 4 ,l/ = (27) 

One is tempted to reduce this equation to a set of ordinary differential equations as before, 
giving 

V{<f>) = ^ 2 (t)s 2 (t)0 2 + \\(t)G 4 m\ (28) 

where the coefficients are running couplings obeying first order differential equations as in 
the massless case. 

However, this is not correct. By considering small excursions in field space, one does not, 
as in the massless case, reproduce the unimproved one-loop potential. This is not surprising. 
In the massless theory, the only scale is set by 0, and thus all logarithms must be of the form 
t = ln((j) 2 / M 2 ) . In the massive theory, there is another scale, and there will be logarithms 
of the form ln((— fi 2 + 3A0 2 )/M 2 ). Thus one cannot easily sum all of the leading logarithms. 
In addition, the scale dependence of the constant term in the potential (the cosmological 
constant) can be relevant. 



In earlier work (and in the review of Ref. ||75||), it was argued that the bounds only depend 



on the structure of the potential at large <ft, and thus the mass term and constant term are 
irrelevant. However, in going from A to the Higgs mass, the structure of the potential near 
its minimum is important, and thus using the naive expression above is not as accurate 
(although it is fairly close). This will be discussed more in the next section. 



More recently, Bando, et al. |77| and Ford, et al. |[7q| , following some earlier work by 
Kastening ||79|| , found a method of including the additional logarithms found in the massive 
theory. In general, they showed that if one considers the L-loop potential, and runs the 
parameters of that potential using L + 1 beta and gamma functions, then all logarithms 
will be summed up to the Lth-to-leading order. The standard model potential, including all 
leading and next-to-leading logarithms, is then (in the 't Hooft gauge) 

V(cf>) = + W + Thil^ 2 (In fr-f) + \Z 2 (In £ - §) 

+ \H 2 (In & - §) + \G 2 (In £ - §) - 3T 2 (in £ - §) ] (29) 
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with W = g 2 (j) 2 /4:, Z = (g 2 + g' 2 )(j) 2 /4:, H = -/i 2 + 3A0 2 , G = -/j, 2 + A0 2 and T = h 2 (f) 2 /2. 
All of the couplings in this potential run with t = hx<p/M. Use of two-loop beta and gamma 
functions will then give a potential in which all leading and next-to-leading logarithms are 
summed over. It was shown by Casas, et al. ]80| that the resulting minima and masses are 
relatively independent of the precise choice of M, as long as this potential is used (use of 
earlier potentials was inaccurate due to a sensitive dependence on the choice of scale). It is 
this potential that will be used to determine bounds on the top quark and Higgs masses in 
the next section. 

First, one should comment on the gauge-dependence of the potential. Bounds on the 
masses of the top quark and Higgs boson are physical quantities, so how can one draw 
conclusions based on a gauge-dependent potential? It has long been known |75| that the 



existence (or non-existence) of minima of the potential are gauge-independent; an early 
calculation of the mass of the Higgs boson in the Coleman- Weinberg model |8T| to two-loops 



in the gauge showed that the gauge-dependence drops out in the final result. A detailed 
analysis of the gauge-dependence of the bounds on the Higgs and top quark masses has 
been carried out by Loinaz, Willey, et al. They find a gauge-invariant procedure 



for determining the bounds, and find that the final result is numerically very close to the 
procedure discussed below. In a model with stronger gauge couplings, however, the gauge 
invariant method might give significantly different results. 

2. Bounds on the top quark and Higgs masses 

The first paper to notice that fermionic one- loop corrections could destabilize the effective 
potential was by Krive and Linde [jS5|| , working in the context of the linear sigma model. 
Later, independent investigations by Krasnikov ||86|| , Hung ||87|| , Politzer and Wolfram |88| 



and Anselm p9[ all looked at the one-loop, non-renormalization group improved potential 



of Eqs. [L9 and 21., and required that the standard model vacuum be stable for all values 



of d>. The first of these was that of Krasnikov ISHI who noted that the bound would be of 
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0(100) GeV, rising to 0(1000) GeV if scalar loops were included. The works of Politzer 
and Wolfram |88| and Anselm |89| gave much more precise numerical results, but ignored 



scalar loop contributions-thus they obtained upper bounds of 80 — 90 GeV on the top quark 
mass. Hung |87[] gave detailed numerical results and did include scalar loops, thus his upper 
bound ranged from 80 GeV to 400 GeV as the Higgs mass ranged from to 700 GeV. 

All of these results are unreliable because the potential used is not valid for large values of 
0. In these papers, the instability would occur for large values of 0, and thus ln(0/cr) is large 
enough that only a renormalization group improved potential is reliable. The first attempt 
to use an improved potential was the work of Cabibbo, Maiani, Parisi and Petronzio fl90 



They included the scale dependence of the Yukawa and gauge couplings, and required that 
the effective scalar coupling be positive between the weak scale and the unification scale. 
Although they didn't use the language of effective potentials, this procedure turns out to 
be very close to that used by considering the full renormalization group improved effective 
potential. Similar results, using the language of effective potentials, was later obtained by 
Flores and Sher [PHI . 



Use of the renormalizat ion-group improved potential will weaken the bounds. The beta 
function for the top quark Yukawa coupling is negative, and thus the coupling falls as the 
scale increases. Thus, the effects of fermionic corrections will decrease at larger scales. 
Compared with the bounds that one would obtain by ignoring the renormalization-group 
improvement, the decrease in the Yukawa coupling at large scales will weaken the upper 
bounds. This effect is not small; the Yukawa coupling for a quark will fall by roughly 
a factor of three between the weak and unification scales. Note that for additional lep- 
tons, the Yukawa coupling does not fall significantly, thus the bounds obtained by the 
non-renormalization-improved potential will not be greatly changed. 

The first attempt to bound fermion masses using the full renormalization group improved 
effective potential (earlier works, for example, never mentioned anomalous dimensions) was 
the 1985 work of Duncan et al. IB3I. Their results, however, used tree level values for the 



Higgs and top masses, in terms of the scalar self-coupling and the "MS* Yukawa coupling" , 
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and found a bound which, to within a couple of GeV, can be fit by the line 

M top < 80 GeV + 0.54:M Higgs (30) 

As we will see shortly, however, corrections to the top quark mass can be sizeable, as much 
as 10 GeV. A much more detailed analysis, using two loop beta functions and one-loop 
corrections to the Higgs and top quark masses (defined as the poles of the propagator), 
was carried out in 1989 by Lindner, Zaglauer and Sher |93], and followed up with more 



precise inputs in 1993 by Sher ||94|| . In all of these papers, the allowed region in the Higgs- 
top mass plane was given-the allowed region was always an upper bound on the top mass 
for a given Higgs mass, or a lower bound on the Higgs mass for a given top mass. The 
allowed region depended on the cutoff A at which the instability occurs. For example, if 
the instability occurs for values of above 10 10 GeV, then one concludes that the standard 
model vacuum is unstable IFF the standard model is valid up to 10 10 GeV (should the 
lifetime of the metastable vacuum be less than the age of the Universe, one would conclude 
that the standard model cannot be valid up to 10 10 GeV). Thus, all of the bounds depend 
on the value of A. 

In the above papers, the effective potential used was the renormalization group improved 
tree-level potential, Eq.|28|. As discussed in the previous section, this would be as precise 
as the precision of the beta functions and anomalous dimensions (two-loop were used) if 
the only logarithms were of the form ln(-^2-); the resulting potential is exact in terms of 
the beta and gamma functions. However, when scalar loops are included, terms of the 
form ln( M ) arise, and these terms are not summed over. In the earlier papers, it was 
argued that when is large, the scalar terms are effectively of the form ln(^j-), and thus the 
difference is irrelevant. But, in determining the Higgs boson mass in terms of the potential, 
the structure of the potential at the electroweak scale is relevant, and thus the difference in 
the form of the scalar loops is relevant. It turns out that this difference is especially crucial 
when the value of A is relatively small (1 — 10 TeV), and less important when A is large 
(10 15-19 GeV), thus the results of the above papers are valid in the large A case. 
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To include the proper form of the scalar loops, one must use the form of Ford, et al. 
discussed in the last section. This analysis was carried out very recently by Casas, 
Espinosa and Quiros |95| and by Espinosa and Quiros ||96|| p|. A very pedagogical review of 



the analysis can be found in Espinosa's Summer School Lectures ||99|| . We now briefly review 
that analysis and present their results. 

Consider the tree level renormalization group improved potential, Eq. At large values 
of (p, the quadratic term becomes negligible, and the question of whether the standard model 
vacuum is stable is essentially identical to the question of whether X(t) ever goes negative. 
If X(t) goes negative at some scale A, then the instability will occur at that scalef]. 

Casas, et al. [p0| , p5| analyzed the question using the full one-loop renormalization group 
improved potential, with two-loop beta and gamma functions, of Eq. They showed that 
the instabilty sets in when A becomes negative, where A is slightly different from A: 

* - * - tL^ (- H - ^ ( k 9,2,2 ( in ^ - sH 

All that remains is to relate the parameters in the potential to the physical masses of the 
Higgs boson and of the top quark. 

It is not a trivial matter to extract the Higgs and top quark masses from the values of 
h(t) and A(t) used in the potential. One can write 

m to p(» = mtop e (l + <WO)) = T/^rW 

m E (ji) =<(HfeM) = ^) (32) 
where the pole masses are the physical masses of the top and Higgs, and <5 top (/i) is the 



radiative corrections to the MS top quark mass. Note that the physical Higgs mass is 
NOT simply the second derivative of the effective potential, since the potential is defined at 



2 See Altarelli and Isidori [98 for a similar and independent analysis. 

3 For a discussion of the relationship between the location of the instability and the required onset 



of new physics, see Ref. ]HHJ. 
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zero external momentum and the pole of the propagator is on-shell; 8g(fj,) accounts for the 
correction. 

The correction <5t op (/z) receives contributions from QCD, QED and weak radiative effects, 
with the QCD corrections being the largest. The QCD corrections have been calculated to 
0(g%) in Ref. | |101| and to 0(gf) in Refs. |102| , |103| , the other corrections were determined 



in Refs. [|104H106| . The correction can be found in Refs. |80|J107|| . The detailed 

expressions for these quantities, which correct several typographical errors in the published 
works, are summarized in an extensive review article by Schrempp and Wimmer [|108|| . The 
largest correction is to the top quark mass; the leading order term is which is 5%, or 
almost 10 GeV. 
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FIG. 2. Perturbativity and stability bounds on the SM Higgs boson. A denotes the energy 
scale where the particles become strongly interacting. 



All of these corrections were included in Refs. [95] and [96|, and reviewed in Ref. |97|. If 
one requires stability of the vacuum up to a scale A, then there is an excluded region in the 
Higgs mass-top mass plane. The result, for various values of A, is given in Figure 2. This 
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figure, in addition, also includes the region excluded by the requirement that the scalar and 
Yukawa couplings remain perturbative by the scale A; these bounds will be discussed in the 
next section. The lower part of each curve is the vacuum stability bound; the upper part 
is the perturbation theory bound. The excluded region is outside the solid lines. Thus, for 
a top quark mass of 170 GeV, we see that a discovery of a Higgs boson with a mass of 90 
GeV would imply that the standard model vacuum is unstable at a scale of 10 5 GeV, i.e. if 
we live in a stable vacuum, the standard model must break down at a scale below 10 5 GeV. 
The curves in Figure 2 are approximately straight lines in the vicinity of M top ~ 170 GeV, 
thus the top mass dependence can be given analytically f97|. For A = 10 19 GeV, we must 
have 

M H (GeV) > 133 + 1.92(M top {GeV) - 175) - ^s " 3 ^" 12 (33) 
and for A = 1 TeV, 

M H (GeV) > 52 + 0M(M top {GeV) - 175) - 0.5 "^p" (34) 



It is estimated |9l)| , |9~T| that the error in the result, primarily due to the two- loop correction 
in the top quark pole mass and the effective potential, is less than 5 GeV. In Figure 3, 
the stability and perturbation theory bounds are given explicitly as a function of A for 
M top = 175 GeV. 
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FIG. 3. Perturbativity and stability bounds on the SM Higgs boson as a function of A for 
M top = 175 GeV. 

Of course, it is not formally necessary that we live in a stable vacuum^. Should another 
deeper vacuum exist, it is only necessary that the Universe goes into our metastable vacuum 
and then stay there for at least 10 billion years. A detailed discussion of the finite tem- 
perature effective potential and tunnelling probabilities is beyond the scope of this review; 



the reader is referred to Refs. |95j and |97[ for the details, as well as a comprehensive list 
of references. In short, the bound in the above paragraph for A = 10 19 GeV weakens by 8 
GeV, and for A = 1 TeV, weakens by about 25 GeV. In all cases, the bound obtained by 
requiring that our vacuum have a lifetime in excess of 10 billion years is weaker than the 
bound obtained by requiring that the Universe arrive in our metastable vacuum. 

We now turn to the question of vacuum stability for models with additional quarks and 



This was first pointed out in Ref. [109 1 
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leptons. 



3. Vacuum Stability Bounds and Four Generations 

In this review, we are concentrating on two possibilities for quarks and leptons beyond 
the third generation: the chiral case and the vector-like case. In the latter case, the quarks 
and leptons cannot couple to Higgs doublets. They will thus have no effect on the vacuum 
stability bounds (except very weakly through their effects on the two-loop beta functions of 
the gauge couplings). We conclude that there are no vacuum stability bounds on the masses 
of vector-like quarks and leptons in models with Higgs doublets. Should one include a Higgs 
singlet, of course, then the vector-like quarks and leptons would couple, and a bound could 
be found on their masses which depends on the singlet Higgs mass as well as the fraction 
of the mass which comes from the singlet vev (since a bare mass is possible). There is one 
exception to this conclusion. As noted by Zheng |110|| , if one adds a vectorlike doublet and 



one or more vectorlike singlets, then Yukawa couplings can exist. He studied the stability 
bounds in that case (using the tree-level potential and one-loop beta functions), assuming 
that the Yukawa couplings were unity, and found that the bounds are much more stringent 
than in the three generation case — the lower bound from vacuum stability and the upper 
bound from perturbation theory come together at a scale well below the unification scale. 

In the chiral case, more specific bounds can be found. It is clear from Eq. that one 
can naively just replace the h 4 term with a summation over all quark and lepton Yukawa 
couplings. This amounts to replacing M top with 

(A4 + M 4 d + l -M 4 E + l -M^ V4 (35) 

where U, D, E and N are additional [/-type quarks, D-type quarks, charged leptons and 
neutrinos, respectively. The | factor is a color factor. This replacement was noted by Babu 
and Ma ||111|| who simply rewrote Eq. by substituting M top with the above expression 
(their paper was written when the top quark mass was believed to be 40 GeV, so the top 
mass was not included in the above.) 
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As discussed above, however, this procedure is not particularly accurate. A more detailed 



analysis, using one-loop beta functions, was carried out by Nielsen et al. ||1 1 2|| For simplicity, 
they assumed that the fourth generation fermions all had a common mass, M 4 ; since the 
quarks must have very similar mass, relaxing this assumption will not significantly affect 
the results. If one assumes that the standard model is valid up to A ~ 10 15 GeV, then they 
found that there is an upper bound on M4 of only 110 GeV. It is easy to see why this bound 
is so stringent. Suppose that M 4 were equal to M top . Then the expression in the above 
paragraph (ignoring the leptons) would be 2M top , and thus the quark contribution would be 
3M top . One might expect the bounds to be smaller by a factor of roughly 3. In fact, it isn't 
quite that severe since the upper line in the allowed region of Figure 2 is not significantly 
affected by the presence of additional generations. Nonetheless, the bound is quite stringent. 



Nielsen et al. ||112|| argued that CDF bounds ||113| , |114|| on stable, color triplets, as well 



as results from the successful description of top quark decays (which occur at the vertex), 
rule out M 4 up to 139 GeV, and thus the standard model cannot be valid up to 10 15 GeV. 
They found that new physics had to start below approximately 10 10 GeV. However, this 
argument has a flaw. It is true that the CDF bounds on stable, color triplets rule out heavy 
quarks with decay lengths greater than a meter or so, and that the successful description of 
top quark decays rule out heavy quarks which decay at the vertex (at least up to about 150 
GeV), but there is still a window for decay lengths in which the quarks would have evaded 
detection. Depending on mixing angles, these decay lengths are quite plausible. Nonetheless, 
theirs is the most detailed analysis to date of the case in which M4 is above 140 GeV (and 
thus the scale of new physics is well below the unification scale). 

A much more detailed analysis was carried out by Dooling et al. ||115| , |116|| They per- 



formed a complete analysis, using the full, two-loop analysis of Casas, Espinosa and Quiros 
(discussed in the last section). They only consider the case in which the standard model is 
valid up to the unification scale, and thus only look at the case in which M4 is very light, 
typically less than 120 GeV. Their work is thus complimentary to that of Nielsen et al. They 
find that the point where the triviality bound and vacuum stability bound come together 
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(see Figure 3) is (for A ~ 10 19 GeV) M 4 < 110 GeV. Thus, if the standard model is valid up 
to the unification scale, only a narrow region of masses still exists between the LEP lower 
bound (roughly 1/2 the center-of-mass energy) and the vacuum stability bound of 110 GeV. 

These works did assume that the quark and lepton masses were all degenerate with a 
mass M 4 . If one relaxes this assumption, then one approximately can replace (8/3) 1 / 4 M 4 
with Eq. [35[ Clearly it is easier to accommodate heavy leptons and neutrinos than heavy 
quarks. 

Hung and Isidori [|117|1 relaxed the assumption of a common M 4 and simply assumed a 
doublet of degenerate quartks with mass Mq and a doublet of degenerate leptons with mass 
Ml. They found that with M L ~ M w , Mq can be extended to 150 GeV before a Landau 
pole appears at the Planck mass. As Mi is raised, Mq should correspondingly decrease if 
one requires that the Landau pole appear at or above the Planck mass. 
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C. Perturbative Gauge Unification 



The use of the RG equations as a tool to set bounds on and, in particlular, to "predict" 
particle masses is an "old" subject. The discussion on the bounds has been carried out in 
previous subsections. This subsection concentrates instead on the use of the RG equations 
to "predict" various masses. In particular, we shall pay special attention to the masses of 
any extra family of quarks and leptons. This analysis only works for chiral fermions. Vector- 
like fermions, having no Yukawa coupling to the SM Higgs field, will not have the desired 
influence on the evolution of the couplings as we shall see below. 

In order to use the RG equations to make "predictions" on the masses, one has to invoke 
either some experimental necessities or some theoretical expectations -or rather prejudices- 
such as fixed points, gauge unification, etc. We shall describe below these concepts along 
with their consequences. To begin, we shall list the RG equations at two loops for the 



minimal SM with three generations |123| 



16vr 2 ^ = 24A 2 + AX(3g 2 - 2.25g 2 2 - 0A5gf) 
-6gf + (16vr 2 )- 1 {30^ 6 
-I3gt + 2gt-80g 2 3 g?]\ 

-lU\ 2 g 2 - 312A 3 - 32g 2 gf} (36a) 

da 2 

(87r 2 )- 1 [-12^ + 6A 2 + ^ 2 

(-12A + 362 2 )-10% 3 4 ]} (36b) 

da 2 

16vr 2 ^ = ^{(41/5) + (167r 2 )- 1 [(199/25)^ 1 2 + (27/5)^ 2 + 

(88/5)s 3 2 -(17/5k 2 ]} (36c) 

16tt 2 M = s 2 4 {-(19/3) + (16tt 2 )- 1 [(9/5K 2 + (35/3)^ 2 + 

24^ -3(? 2 ]} (36d) 
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16vr 2 ^ = s 3 4 {-14 + (^-MCH/S)^ + % 2 

-52^3 2 -4<? 2 ]} (36e) 

To set the notations straight, our definition of gi uses the SU(5) convention and is related 



to the U(l)y coupling g' by g\ = J 3/5</. In the evolution of these couplings we will neglect 
the contributions coming from the lighter fermions. 

In the above RG equations, clearly the important couplings are those of the top quark 
Yukawa and of the Higgs quartic couplings, and, to a certain extent, also the QCD coupling. 
As we have seen earlier, one important use of such equations is by following the evolution of 
A with the initial value of gt fixed by its experimental value. Requiring A to be positive (for 
vacuum stability reason) at least up to the Planck scale allows us to set a lower limit on the 
Higgs mass to be ~ 136 GeV. This use of the RG equations is rather solid in the sense that 
it relies only on the vacuum stability criterion of quantum field theory. Other uses which 
are discussed below are more speculative but are quite interesting in that several predictions 
can be made and can be tested. 

In dealing with RG equations, a natural question that comes to mind is whether or not 
there exist stable fixed points. Basically, a stable fixed point is a point in coupling space to 
which various couplings converge regardless of their initial values. This is an attractive idea 
that has wide applications in many fields of physics, such as critical phenomena- to mention 
just one of many. In particle physics, there were many speculations concerning the nature 
of such fixed points if they truly exist. For example, Gell-Mann and Low | 124|| , and sub- 



sequently, Johnson, Wiley and Baker ||125|| have speculated that quantum electrodynamics 
might possess an ultraviolet stable fixed point which will render QED finite. Other more 
"recent" speculations dealt with the very interesting subject of the origin of particle masses- 
at least of the heavy one(s). 

In general, a stable fixed point appears as a zero of the (3 function which would be 
meaningful only if one has a full knowledge of such a function. In the absence of such a 
knowledge, one might have to resort to approximations allowed by perturbation theory. In 
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regions where various couplings can be considered to be "small" enough so that the use of 
one or two-loop (3 functions might be justified, one would try to "run" the couplings over a 
large region of energy and see if they converge to a point for arbitrary initial values. If such 
a point is found, say by a numerical study of the RG equations, one would qualify this as 
a fixed point. Such an approach has been pioneered by Pendlenton and Ross |126f , and, in 



particular by Hill |127|| , where the fixed points are of the infrared nature. Of relevance to 



this report is the suggestion by various authors that a fourth generation might be needed 
for the existence of such a fixed point. 

Let us first summarize what has been done for the top quark mass and subsequently 
describe works related to the masses of a fourth generation. 

Pendleton and Ross ||126|| were the first to suggest a relationship between the top quark 
Yukawa coupling and the QCD coupling as a result of an infrared (IR) stable fixed point. 
To see this, one can combine the one-loop RG equations for g t and the QCD coupling 
(first terms on the right-hand side of Eq. (38)) to form a RG equation for the ratio g t j gz, 
namely 

1C 2^(9t/gs) 9 2 2 3. 2 , /2\ 2 /2 . . 

— dt — = 2 9t ~ 9s ~ V 9 9 >~3 9 ' ( > 

Ignoring the electroweak contributions in Eq. (39), there is an IR fixed point obtained by 
setting the right-hand side to zero. Pendleton and Ross obtained a relation 

9t = 9t, irs = \gl (38) 

It turns out that the above relation gives too low of a mass for the top quark. In fact, the 
original prediction ||126|| using Eq. (40) and a value of a 3 = gf/47r ~ 1/7 (at a scale of 



~ 2M\y) gives a mass of ~ 110 GeV. Using the current value of 0:3 ~ 0.12 (at the Z mass), 
the prediction would have been even lower, even after electroweak corrections are included. 
It goes without saying that this cannot be true for we already know that the top quark 
mass is ~ 175 GeV. As pointed out by Hill | |127| | long before the discovery of the top quark, 



the Pendleton-Ross fixed point is only a quasi fixed point in the sense that it can never be 
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reached at the scale of interest ~ 100 GeV. Hill proposed an intermediate fixed point that 
can be found by setting the f3 function in the RG for g t to zero with a "slowly varying" g 3 
replaced by a constant taken to be some average value between two scales: 100 GeV and 
10 15 GeV. This translates into an approximate relation 

\g 2 t « m 2 , (39) 

With g 3 2 ~ 0.7, Hill made the prediction for the top quark mass to be ~ 240 GeV. This is 
now known to be much too large, although at the time the prediction was made, it appeared 
to be a plausible value. 

In the above discussion for the top quark mass as a result of an IR stable fixed point, 
one feature clearly emerges: a heavy fermion is needed to drive the evolution toward a 
fixed point. This point was made even clearer in a detailed study of Bagger, Dimopoulos 



and Masso |[L28|| . These authors made two assumptions: the first one is the existence of a 
desert between the weak scale Mw and some Grand Unified scale Mx ~ 10 15 GeV and the 



second one being that of perturbative unification. The question asked in Ref. [12S] was the 



following: what should the initial values of various Yukawa couplings at Mx be in order for 
those couplings to reach the fixed in a "physical time" tw = ^^h^-p^) ~ 1/5? Again, it 
turns out that "large" initial values (at Mx) of Yukawa couplings guarantee that the fixed 
point is reached in "physical time". What is this fixed point and what does it say about 
masses of possible extra generations if they exist? We shall describe below the salient points 
of the analysis of Ref. [|1 28| . 



We begin the discussion of Ref. J128| ] with the following one-loop RG equations for the 
quark and lepton Yukawa couplings: 

dT c 



dt 
dT L 



2(G Q -T)T Q -3Tr(S 2 u ), (40a) 
2(G L -T)T L -3Tr(S 2 E ), (40b) 



dt 

where various Yukawa factors are defined as Ty = Tr(Y'Y) with Y = U, D, E or N, Tq = 
T v + T D , T L = T E + T N , T = 3T Q + T L , S v = WU - DW and S E = E^E - N*N. Finally 
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the gauge factors G's are defined as Gjj = Gp, = Gq = 8g>| + \g\ and Ge = G^ = Gl = \g\ 
where the contribution from g\ has been neglected. Notice that t = ln( ; 
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FIG. 4. (a) The trace Tq at M\y as a function of the trace Tq and Mx for Np = 8 and Tj_, = 0. 
The dotted line denotes the radial quark fixed point. For Tq(Mx) > 0.1, the fixed point is reached 
in physical time. (b)Ti(Mv^) as a function of Tl(Mx) for Np = 8 and Tq = 0. 



To simplify the discussion, Ref. |[L28|| first assumed degenerate quarks and degenerate 
leptons so that Sy = as well as Sp = 0. If the gauge couplings can be approximated as 
constant (or very slowly varying), Gq and Gl in Eqs. (42) can be replaced by some averages 
similar to the procedure used by Hill ||127| . Let us denote these averages by Gq and Gl- 
It is then easy to see that Eqs. (42) have the following two distinct fixed points: Gq = T 
(quark radial fixed point) and Gl = T (lepton radial fixed point). Whether or not these 
fixed points are reached will depend on the initial values of Tq or Tl at the Grand Unified 
scale M X - If T{M X ) is below some critical value T(M x ) m i n , the fixed point G will not be 
reached in "physical time" t ~ 1/5: the value of T at My/ will be less than the fixed point 
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value G = T. This fact has allowed Ref. ||128|| to set an upper limit on heavy fermion masses 



with the upper limit being the IR stable fixed point. Setting Sy = T L = 0, Ref. ||12S|| plotted 
Tq at the weak scale as a function of its value at Mx- This is shown in Fig. 4a. The result 
for the lepton case is shown in Fig. 4b. 

The figure shows the result for eight families. We are, of course, concerned only with 
four families which are still allowed, within error, by precision electroweak results. For four 



families, Ref. [|128|1 gave the following upper bounds on T: 

T Q < 2.7, (41a) 



T L < 3.4. (41b) 

(The above numbers used values of gauge couplings which are now outdated.) This translates 
into the following distinct bounds on fermion masses for four families: 

EM| < (290 GeK) 2 , (42a) 

EM| < (325 GeV) 2 . (42b) 

(These bounds would be slightly less if recent values of the gauge couplings are used.) The 
above bound for the quarks, for example, would translate roughly into a bound on the mass 
of a degenerate fourth generation quark (after subtracting out the top quark) as M UtD < 
164 GeV. Is this what one should be aiming for when one tries to look for fourth-generation 
quarks? 
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FIG. 5. Flow of A and gt towards fixed points in the standard one-Higgs doublet model. Open 
circles denote initial points. Crosses denote final fixed points. 



Hill, Leung and Rao [|1 29j] made an extensive study of the RG fixed points and their 
connections with the mass of the Higgs boson(s) for the one-Higgs doublet and the two-Higgs 
doublet cases, and for up to five generations. In this work, the Higgs quartic coupling(s) 
is run simultaneously with the various Yukawa couplings and, as a result, one clearly sees 
again the interplay between heavy fermions and the Higgs field. This is shown for example 
in Figs. 1 and 2 of Ref. ||129|| for the one-Higgs doublet case with three generations, which 
are reproduced in Figs. 5 and 6. 
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FIG. 6. Relation between the Higgs mass and the top quark mass in the standard one-Higgs 
model. 

Other attempts of using the RG fixed points to construct fermion mass matrices have 



been made, for example in Ref. | 153 |. However the "predicted" value for the top quark mass 
is now outdated. 



A different approach was taken by one of us (P.Q.H.) concerning the influence of a 
possible fourth generation on the evolution of all couplings of the SM and not merely the 
Yukawa couplings. In particular, the question that was asked was whether or not there can 
be gauge unification in the nonsupersymmetric SM and under which conditions this can be 
achieved. As we shall see below, it turns out that a fourth family of chiral fermions will be 
needed and that their masses are found to be fairly constrained. 

The possibility of coupling-constant unification of the three gauge interactions of the 
Standard Model (SM) is, without any doubt, one of the most important issues in particle 
physics. Coupling-constant unification is a necessary, but not sufficient, condition for a 
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Grand Unification of the SM 1 11 3 1| 1 1 3 3|] . Such a possibility is particularly attractive since 



it would provide a unified explanation for a number of puzzling features of the SM such as 
electric charge quantization for example. 

There are various ways that the three gauge couplings can get unified. The simplest way 
is to assume that there is a "desert" (i.e. no new physics) between the electroweak scale and 
the scale at which unification occurs. Simply speaking, the three gauge couplings are left to 
evolve beyond the electroweak scale under the assumption that there is no additional gauge 
interactions of the type which would modify the evolution of some of the gauge couplings. 
A more complicated way would be to assume that there is one or several intermediate scales 
where partial unification among two of the three couplings occurs. This might well be 
the case. However we shall restrict ourselves, in this report, to the simplest scenario of 
unification with a "desert" and search for conditions under which this can be achieved. This 
was the approach taken by one of us (P. Q. H.) [|6^] . 

We will proceed in two steps. First, we will present the evolution of the gauge couplings 
and show the places where they cross, ignoring any heavy threshold effects that might be- 
and should be if there truly is unification- present. In this discussion, we will show both the 
minimal SM with three generations and the one with an extra fourth generation. We shall 
see that, under certain restrictions on the masses (fourth generation and Higgs masses), 
the latter possibility provides a better "convergence" of the three gauge couplings. By 
convergence under the quotation marks, we mean that they do not precisely meet at the 
same point. The "true" convergence will be shown to be accomplished by the inclusion of 
heavy threshold effects. In fact it would be senseless to claim unification without taking into 
account such effects. 

We first summarize the situation with three families. 

The first task is to integrate Eq.(38) numerically and look for the places where the 
couplings meet, disregarding for the moment the possibility that there might be unification. 
We then have to set up some kind of criteria to decide on how close to each other all 
three couplings have to be in order for them to have a chance of actually converging to 
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a single scale, once heavy particle thresholds, such as those of the X and Y bosons of 
577(5) for instance, are taken into account. Once these criteria are satisfied and heavy 
particle threshold effects are included, one can put an error on the unification scale and, 
consequently, an error on the proton lifetime. 

Fig. 7 shows the evolution, without heavy threshold effects, of g 3 , g 2 , and gi of SU(3) <g) 
577(2) <S> U(l) for the case with three generations. Clearly these three couplings do not 
converge. 

1.50 i i i i i i i i i i i i i i i i i 
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FIG. 7. Evolution of couplings in the three generation case. 
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The question is: How far apart are they from each other and at what scales? As far as 
the scales are concerned, we will be interested only in those which are above some minimum 
value implied by the lower bound on proton decay. A rough estimate of that lower bound is 
obtained by noticing that T p ^ e+7T o(yr) m 10 31 (M G /4.6 x 10 14 ) 4 . This gives M G > 1.3 x 10 15 
GeV (corresponding to ln(£'/175) = 29.64 on the graph) for T p ^ e +^o{yr) ^ 5.5 x 10 32 . The 
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next question is the following: Starting from M™ n ~ 1.3 x 10 15 GeV, how far apart are the 
three gauge couplings from each other at a given energy scale? As stated above, the reason 
for asking such a question stems from the fact that, if the SM were to be embedded at Mg in 
a Grand Unified model such as 577(5) [|132j| (for instance), the decoupling of various heavy 
GUT particles would shift the three couplings from a common «g to (possibly) different 
values. As shown below, for a wide range of "reasonable" heavy particle masses, such an 
effect produces no more than ~ 5% shift from the common value and in the same direction. 
It turns out that the modified couplings can differ by no more than ~ 4%. From this 
a reasonable criterion would be to require that, at a scale Mq > rriQ tn , the three gauge 
couplings are within 4% of each other. 

We shall take 577(5) [|132j| as a prototype of a Grand Unified Theory. Let us assume 
the following heavy particle spectrum: (X, Y) = (3, 2, 5/6) + c.c. with mass My, real scalars 
(8, 1, 0) + (1, 3, 1) + (1, 1, 0) (belonging to the 24-dimensional Higgs field) with mass M 2 4, and 
the complex scalars (3, 1, —1/3) (belonging to the 5-dimensional Higgs field), with mass M 5 . 
(The quantum numbers are with respect to SU(3) Cg> SU(2) ® U(l).) The heavy threshold 
corrections are then 111 3411: 



A ' = ?'"(§ i '-^r" 1 (ir) + A ™ - < 43a > 

47T My oOlT M 5 

^-1 + >^-^wJ +a " ro < (43b) 



where 



with ki = 1/2,3/2,-1 for i = 1,2,3, is the correction coming from possible dimension 5 
operators present between Mq and Mpi anc k. The modified gauge couplings can be expressed 
in terms of the unified coupling (Xq (at M G ) as: 
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where i = 1,2, 3. We then define the fractional difference between the modified gauge 
couplings as: 



da = ^-A (46) 

for i,j = 1, 2, 3 and the definition refers to oti as being the larger of the two couplings. For a 
wide range of heavy particle masses (in relation with Mq) and the parameter rj appearing in 
Aj, and for ~ 0.024 — 0.028, it is straightforward to see that dij can be at most 4% p3 . 
^From this simple analysis, one can reasonably set a criterion for a given scenario to have a 
chance of having gauge coupling unification: the fractional difference among the three gauge 
couplings at some scale Mq > M™ n should not exceed 4%. 

For the SM with three generations and taking into account the presence of M™ n ~ 
1.3 x 10 15 GeV, one finds the following trend: d^2 decreases from 3% as one increases the 
energy scale beyond M™ n , while d^% increases from 4% and d 2 \ also increases from 7%. 
(For example, at M G ~ 3.3 x 10 15 GeV, d 32 ~ 1.4%, d 31 ~ 8.4% and d 21 ~ 9.7%. /From 
these considerations- and not from just "eyeballing" the curves- one might conclude that 
the minimal SM with three generations does indeed have some problem with unification of 
the gauge couplings. 

There is a drastic change to the whole scenario when one postulates the existence of a 
fourth generation of quarks and leptons [65[. The main reason is the fact that the Yukawa 



contributions to the running of the gauge couplings appear at two loops. In the three 
generation case, the top Yukawa coupling actually decreases sligthly with energy because 
its initial value is partially cancelled by the QCD contribution (at one loop). As a result, 
the presence of a heavy top quark is insignificant in the evolution of the gauge couplings 
at high energies when there are only three generations. The presence of more than three 
generations drastically modifies the evolution of the Yukawa, Higgs quartic self-coupling, 
and the three gauge couplings. For example, with a fourth generation which is sufficiently 



42 



heavy, all Yukawa couplings grow with energy, significantly affecting the evolution of the 
gauge couplings. It turns out, as we shall see below, that the Yukawa couplings can develop 
Landau poles below the Planck scale. If there were any possibility of gauge unification, one 
would like to ensure that it occurs in an energy region where perturbation theory is still 
valid. Furthermore, the unification scale will have to be greater than Mq iu (as discussed 
above). As we shall see, the validity of perturbation theory plus the lower bound on the 
proton lifetime put a severe constraint on the masses of the fourth generation. 

The two-loop renormalization group equations applicable to four generations are given 



by f!23lj65j: 




24A 2 + 4\(3g 2 + 6g 2 q + 2g 2 - 2.2hg 2 2 - 0.45^) 
-2(3^ 4 + 64 + 2^) + (167r 2 )- 1 {30^ 



+48^ + 16gf - [3gf + 6g 4 q + 2gf - 80g 2 3 (g 2 
+2^ 2 )]A - 6A 2 (24^ 2 + 48g 2 q + 16gf) - 312A 3 
-32flg(rf + 2$} 



(47a) 




g 2 {§g 2 + I2g 2 q + Ag 2 - 16^ 2 - A.hg 2 2 - 1.7 g* + 
(Stt 2 )- 1 ^ 4 - 2.25g 2 (6g 2 + 3g 2 + 2gf) 
-12g*-(27/4)gt-3gt + Q\ 2 + g 2 
(-12A + 36£ 2 )-(892/9)s 3 4 ]} 



(47b) 




g 2 {6g 2 + 12g 2 + Agf - 16g 2 - 4.5^ - 1.7 gl + 
(8n 2 )- 1 [3g 4 q -g 2 (6g 2 q + 3g 2 + 2g 2 ) 
-12g\ - (27/4)s 4 - 3gf + 6A 2 + g 2 
(-16A + 40( 7 2 )-(892/9)( ? 4 ]} 



(47c) 



16tt 2 J- = g?{6g 2 + 12g 2 + Agf - A.5(g 2 2 + gf) + 

(Stt 2 )- 1 ^ 4 - g 2 q (6g 2 q + 3g 2 + 2gf) - 12g* 
-(27/4)^ 4 -3^ 4 + 6A 2 -16A^ 2 ]} 



(47d) 
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16n2 it = ^ 163 / 15 ) + (16tt 2 )- 1 [(787/75)^ + 6.6g 2 2 + 
(352/l% 3 2 - 3.4<? 2 - AAg 2 q - 3.6«? 2 ]} 



(47e) 




^ 2 4 {-(ll/3) + (IGtt 2 )- 1 ^ 2 + (133/3),? 2 + 



32« ?3 2 -3< ? 2 -3« ? 2 -2« 7 f]} 



(47f) 




^{-(34/3) + (16tt 2 )- 1 [(44/15)^ + 12^ 2 



(4/3)^-4^-8^]}. 



(47g) 



For simplicity, we have made the following assumptions: a Dirac mass for the fourth neutrino 
and the quarks and leptons of the fourth generation are degenerate SU(2)l doublets. The 
respective Yukawa couplings are denoted by g q and g x respectively. Also, in the evolution 
of the quartic coupling A and the Yukawa couplings, we will neglect the contributions of r 
and bottom Yukawa couplings, as well as the electroweak gauge couplings, g\ and g 2 , to the 
two-loop (5 functions since they are not important. Also, as long as the mixing between the 
fourth generation and the other three is small, one can neglect such a mixing. 

In the numerical analysis given below we shall fix the mass of the top quark to be 175 
GeV. We shall furthermore restrict the range of masses of the fourth generation so that the 
Landau poles lie comfortably above 10 15 GeV, in such a way that unification occurs at a 
scale which would guarantee the validity of perturbation theory as well as satisfying the 
lower bound on the proton lifetime. Concerning the former requirement, it basically says 
that one should look at unification scales where the values of the Higgs quartic and Yukawa 
coulings are still sufficiently perturbative that one can neglect contibutions coming from 
three-loop (and higher) terms to the f3 functions. 

Fig. 8 shows g\, g\ and g\ as a function of energy for a particular set of masses: itlq = 151 
GeV, rriL = 95.3 GeV, where tuq and rriL denote the fourth generation quark and lepton 
masses respectively. 
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FIG. 8. Couplings as a function of energy for the set of masses given in the text. 



It is already well known, from the discussion in the previous section, that, by adding 
more heavy fermions, the vacuum will tend to be destabilized unless the Higgs mass is large 
enough. As we have seen, the vacuum stability requirement is equivalent to the restriction 
A > 0. Furthermore, the heavier the Higgs boson is, above a minimum mass that ensures 
vacuum stability, the lower (in energy scale) the Landau pole turns out to be. It turns out 
that this Landau pole should not be too far from M™ n otherwise g^, gi and g\ would not 
come close enough to each other. On the other hand, it should not be too close either because 
of the requirement of the validity of perturbation theory. These considerations combine to 
give a prediction of the Higgs mass, namely m# = 188 GeV for the above values of the 
fourth generation masses p5| . The dependance of the Higgs mass on the fourth generation 
mass in this analysis is obviously striking. 

Following the criteria that we have set for taking into account the heavy threshold effects, 
the midified couplings «j(Mg) expressed in terms of a^Mc) (which can be read off from the 
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graph) and the threshold correction factors A, are given by: I/ol^Mq) = 1/oh(Mg) + A*. 
The choice of the mass scales M 5 , M 24 , My, and the parameter rj is arbitrary and is only 
fixed to a certain extent by the requirement that a^MaYs should be as close to each other 
as the precision allows. As an example, the choice M5 = M G , M 24 = M G , My = 0.5M G 
and 77 = 10 (where we have picked M G pa 3.5 x 10 15 GeV) transforms a^{M G ) = 0.0278, 
a 2 (M G ) = 0.0273 and «i(M G ) = 0.0285 (values that can be read off Fig. (8)) to a 3 (M G ) = 
0.02735, a 2 (M G ) = 0.02662 and ai(M G ) = 0.02705. From these values, one can conclude 
that the couplings are practically the same with all three equal to a G ~ 0.027 or l/a G pa 37. 

The above simple exercise simply shows that, with just an additional fourth generation 
having a quark mass niq pa 151 GeV, a lepton mass mi pa 95.3 GeV and a Higgs mass 
niu ~ 188 GeV, unification of all three gauge couplings in the nonsupersymmetric SM 
can be achieved after one properly takes into account threshold effects from heavy GUT 
particles ||65|| . Other combinations of masses are possible for gauge unification but their 



values will not be much different from the quoted ones, the reason being the requirement 
that the mass range of the fourth generation be restricted to one that will have Landau poles 
only above 10 15 GeV. Do the masses given above satisfy the requirement of perturbation 
theory? In fact, at the unification point M G = 3.5 x 10 15 GeV, one has (with a, = gf/Air): 
at = 0.4, a q = 0..16, a\ = 0.48 and X/4-k = 0.19. Although these values are not "small", 
they nevertheless satisfy the requirements of perturbation theory, namely ott, q ,i ~ 1 and 
A/4-7T < 0.4. (The latter requirement comes from lattice calculations which put an upper 
bound on the Higgs mass of ~ 750 GeV.) For comparison, a t in the three- family SM has a 
value of 0.016 at a comparable scale and this explains why it is unimportnat in the evolution 
of the SM gauge couplings. 

An important consequence of a fourth generation in bringing about gauge unification is 
the value of the unification scale itself. In the example given above, it is M G = 3.5 x 10 15 



GeV ||65|| . In the nonsupersymmetric SU(5) model, the dominant decay mode ofthe proton 
is p — > e + 7T° and the mean partial lifetime is r p ^ e + n o(yr) pa 10 31 (Mg/4.6 X 10 14 ) 4 . Taking 
into account various uncertainties such as heavy threshold effects, hadronic matrix elements, 
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etc., the predicted lifetime is r p _> e + n o(yr) « 3.3 x 10 34±2 to be compared with T^ e+n0 (yr) > 
5.5 x 10 32 |S5j| . Notice that the central value is within reach of the next generation of 
SuperKamiokande proton decay search. 

Another hint on the masses of a fourth generation comes from considerations of models 
of dynamical symmetry breaking a la top-condensate [I6| This will be discussed in Section 
5 where one can see how the original idea of using the top quark as the sole agent for 
electroweak symmetry breaking (in the form of ti condensates) led to a prediction for the top 
quark mass (before its discovery) to be much larger than its experimental value. The original 
form of this attractive idea obviously has to be modified, most likely by the introduction of 
new fermions such as a fourth generation or SU (2)-singlet quarks. 

In the above discussion on perturbative gauge unification, as well as in the subsequent 
related discussion in Section V, the issue of the gauge hierarchy problem is not considered. 
Such an issue is beyond the scope of the largely phenomenological approach that we are 
taking. This point was alluded to in our Introduction where we stressed that none of 
the reasons given for considering quarks and leptons beyond the third generation is fully 
compelling, but each, including the one on perturbative gauge unification, is suggestive. 
It is certainly possible that the "solution" of the gauge hierarchy problem will not affect 
the above arguments; the recently developed alternative to supersymmetry and technicolor, 
TeV-scale gravity, for example, may not appreciably change results on gauge unification. A 
full consideration of the gauge hierarchy problem is beyond the scope of this review. 
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D. Mixing Angles 



In previous sections, we have seen that the masses of quarks and leptons, although arbi- 
trary, are constrained by phenomenological considerations as well as vacuum stability and 
perturbation theory. The mixing angles of quarks and leptons are also arbitrary, however 
there are no constraints from vacuum stability and perturbation theory (and only weak 
phenomenological constraints). Thus, a much wider range of mixing angles can be accom- 
modated, and one can only be guided by considering various models for these angles. In this 
section, we will discuss plausible models for mixing angles. Since we know that the quark 
sector has nonzero mixing angles, but that the lepton sector may not, we will first look at 
the lepton sector, and then the quark sector. 



1. Leptons 

The only phenomenological indication of any mixing in the lepton sector comes from 
neutrino oscillations. At the time of this writing, there are three indications of oscillations: 
solar neutrinos [135], atmospheric neutrinos [136] and LSND [|137|| . It is difficult, although 



not quite impossible, to reconcile all three of these in a three generation model. If there are 
four light neutrinos, in this case, the fourth neutrino must be sterile (an isosinglet) in order 
to avoid the bounds from LEP. Such a neutrino could exist without requiring the existence 
of any additional fermions. It is likely that the situation will be clarified within a year or so 
at Superkamiokande and the Solar Neutrino Observatory. A detailed discussion of neutrino 
oscillations and their phenomenology, including the recent strong evidence for atmospheric 



neutrino oscillations at SuperKamiokande, can be found in Ref. ||138|| . We will defer to that 
review in this paper, and will not discuss the possibility of light isosinglet neutrinos further. 

We certainly will, however, discuss the case in which a fourth generation neutrino is very 
heavy. This will automatically occur if the fourth generation is vectorlike. Even if it is 
chiral, models exist that can give such a mass. Recently, one of us ]|139|| has considered a 
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model of neutrino masses with four generations where one can obtain dynamically one heavy 
fourth generation and three light, quasi-degenerate neutrinos. Ref. ||140|| has also considered 
a scenario with four generations which has similar consequences. 

Suppose that the heavy leptons form a standard chiral family, with a right-handed neu- 
trino. The bounds from the Z width obtained at LEP force the mass of the N and E to be 
greater than 45 GeV. Are there any phenomenological bounds on the mixing? In analogy 
with the quark sector (as well as the prejudice from most models), one expects the mixing to 
be the greatest between the third and fourth generations. This will affect the tv t W vertex, 
multiplying it by cos#, where 9 is the mixing angle. Since all r decays occur through this 
vertex, the result will be a suppression in the overall rates. For some time, it was believed 
that the mass of the r was 1782 ± 2 MeV, and the measured rate was too low; mixing with 
a fourth generation was a straightforward explanation ||141| -143l. However, the r mass has 
now been measured to much higher precision at BES ||144|| to be 1776.96 ± 0.2 ± 0.2 MeV, 
and the measured rate is now in agreement with theoretical expectations. This has been 
analyzed by Swain and Taylor | p.45|jl46|| , who find a model- independent bound on the mixing 
of sin 2 9 < 0.007. A similar bound can be obtained for mixing between the fourth generation 
and the first two, although one expects those angles to be smaller. 

What values of the mixing might one expect? There are four plausible (in the view of 
the authors, of course) values of the mixing angle between the third and fourth generations: 
(a) sin 2 9 = m T /rriE {b) sin 2 9 = m UT /m^ (c) sin 2 9 = m^ /v /m 2 pl (d) sin 2 9 = 

The first of these occurs in typical see-saw models. The second occurs in models in 
which the mixing occurs only in the neutrino mass matrix. The third occurs in models with 
a global or discrete lepton-family symmetry broken by Planck scale effects, and the fourth 
occurs when the symmetry is not broken by Planck scale effects. We now discuss each of 
these. 

The first relation, sin 2 9 = m T /rriE, occurs in models in which the 2x2 mass sub-matrices 
(0 A\ 

are of the form . If the neutrino and charged lepton mass matrices are of that 

\A B 
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form, then the mixing angle is given by ym T /mE — \Jm UT /mN-, which gives sin 2 9 = m T /niE 
for realistic values of the v T mass. Models of this type were pioneered by Weinberg ||147|| and 
Fritzsch [ |148| | , who noticed that they will give the successful relation for the Cabibbo angle: 
sin 2 6 C = md/m s . Fritzsch also showed [14S] that there are some very simple symmetries 
which automatically give this relation. When the rate for leptonic decays was believed to 
be too low, Fritzsch ||1 50|| used this relation to propose that a fourth generation lepton of 
100 — 200 GeV could account for the discrepancy. 

As noted above, there is a lower bound on the mixing between the r and the E given 
by sin 2 6 < 0.007. Using the Fritzsch relation, this becomes a lower bound on m^, which 
is given by wie > 250 GeV. This is very near the bounds from perturbation theory. We 
conclude that a very slight improvement in the uncertainties in the r decay rate will rule 



out the very general relationship sin# = ym r /m£ (or discover the effect!). 

The second relationship, sin 2 ^ = m Ur /rriN, will occur in models in which, because of 
some discrete or global symmetry, the charged lepton mass matrix is diagonal. The Fritzsch 
relationship will then give sin 2 9 = m VT jm^. Given the cosmological bound on the v T mass, 
this gives a value of sin 2 6 which is less than 10 -10 . The E or N lifetime (whichever is the 
lighter) lifetime will then be in the picosecond-nanosecond range, with extremely interesting 
phenomenological consequences. 

Suppose that one simply assumes that a discrete symmetry forbids any mixing at all 
between the E and N and the other three generations. This is simply an extension of the 
familiar electron-number, muon-number and tau-number conservations laws. In this case, 
the mixing angle vanishes and the lighter (the E or the N) is absolutely stable. As will be 
seen in the next Section, this would be cosmologically disastrous if the E is stable, but not 
if the N is stable. 

Finally, one can assume the discrete symmetry which forbids mixing, but note that Planck 
mass effects are expected to violate all discrete and global symmetries. That means that 
higher dimension operators, suppressed by the Planck mass, will violate these symmetries. 
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Two such examples are given by Kossler et al. | 152|| . The mixing angle is then given by 
sin# ~ My/ /Mpi. This gives a lifetime for the lighter of the E or N of approximately ten 
years, which is very near the bound for charged leptons, discussed in the next Section. 



2. Quarks 

In the lepton case, one could obtain stringent bounds on mixing with a fourth generation 
by considering precise measurements of leptonic decays with theoretical expectations. Here, 
such precision (both theoretical and experimental) is impossible. One can still obtain bounds 
on mixing between the first two generations and a fourth from the unitarity of the CKM 



matrix. As noted in the Particle Data Group Tables f40 |, the mixing angle between the 
first and fourth generations, V u r> must be less than 0.08. However, other bounds are much 
weaker ||153| , [L56|| — the mixing angle between the second and fourth generations, V c d, is only 



bounded by sin 2 6 < 0.5. In the top sector, one can use constraints ||153|| from Ki — > 
to find that Re(V* v Vdu) < 8 x 10 -4 . Since these are all mixings between the first two and 
fourth generations, they are expected to be very small-of greater interest is the bound on 
the mixing between the third and fourth generations. The value of the V t b element in the 



CKM matrix is greater than 0.99 (leaving very little room for such mixing), however |157 
this is determined assuming only three generations and CKM unitarity. If one relaxes this 
assumption, the value of V t b could be as small as 0.05 [^0| , |157|| . Thus, the mixing angle 



between the third and fourth generations could be extremely large, and there are effectively 
no phenomenological constraints on such mixing, if the D is sufficiently heavy to avoid 
affecting top quark decays. 

Bounds from K° — K and B° - B mixing are also not very strong. The experimental 
value of B° — B° mixing was the first indication that the top mass might be heavy, and the 
observation that it is, in fact, heavy means that only very weak bounds on fourth generation 
masses and mixings may be obtained. Similarly, the large number of phases and angles in 
the four-generation CKM matrix (3 and 6, respectively), implies that only weak constraints 
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can be found from e'/e. 

The reason that the bounds for these flavor changing processes are so weak is due to 
the GIM mechanism. This mechanism only applies if the fourth generation is chiral. If it is 
composed of vector-like isosinglets or isodoublets, then the GIM mechanism will break down 
and one will have Z-mediated flavor-changing neutral couplings (FCNC). In addition, one 
can have an effect on flavor-diagonal neutral currents (FDNC), since mixing of a doublet 
quark with a singlet will reduce its left-handed coupling. To be more specific, consider the 
case of a Q = —1/3 isosinglet quark, D. This case has been analyzed in great detail by 
Barger, Berger and Phillips ||154|| . The mixing between mass and weak eigenstates is given 
by 
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In the basis where the Q = 2/3 mass matrix is diagonal, the first three rows and column of 
the V matrix are just the usual CKM matrix. The fourth row is not relevant for the weak 
interactions with the W and Z. Since the entire matrix is unitary, the CKM matrix will not 
be, leading to a suppression of flavor-diagonal couplings. The Z couplings are given by 



C-fcnc = \az^Zijq iL ^Z^ L 



(49) 



where, using the unitarity of the 4x4 matrix 

"j 

The FDNC couplings are given by 

C-FDNC = 9zJ2Vi^ Z V 



Zij — Oij - V oi V j 



1 



(50) 



-75) 



- sin 2 9 W 



(51) 



Thus, mixing with the D quarks reduces the direct left-handed FDNC couplings of the light 
quarks, while leaving the right-handed quark couplings unchanged. For &Q = 2/3 isosinglet, 
the results are very similar. 
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Barger, Berger and Phillips (BBP) then analyze a large number of processes, constraining 
elements of the 4x4 matrix. Their paper (written in early 1995) has a huge number of 
references to papers dealing with isosinglet quarks; the reader is referred to it for earlier 
references. Since the top quark was discovered, calculation of the effects of isosinglets became 
much more precise, since a major uncertainty in the calculations had disappeared. This was 
the motivation of BBP for reanalyzing all of the constraints on the elements of the matrix. 
BBP examine Z-decays, meson-antimeson oscillations, Kl, B°, D° —* B,D — > Xl + l~ , 

and radiative B decays, for both the Q = —1/3 and Q = 2/3 cases. They list bounds on 
the off-diagonal matrix elements. 

The above argument that V t b could be very small, while V t o could be very large, does 
not work in the isosinglet Q = —1/3 case. The reason is that the mixing, by reducing 
the left-handed FDNC couplings of the left-handed quarks, can be constrained from high- 
precision SLC and LEP results. In the case where one just adds a chiral family, the FDNC 
couplings are not reduced. This is discussed in detail by Nardi, et al. ||155| 1, who show that 
the precision data gives \V 0d \ 2 < 0.0023, \V 0s \ 2 < 0.0036, \V 0b \ 2 < 0.0020. Unitarity of the 
matrix then gives \V d\ > 0.996, which in terms gives \V q u\ < 0.09 for q = u,c,t. Thus, the 
mixing between the top quark and the D cannot be very large. 

In the Q = —1/3 case, BBP find that the off-diagonal matrix elements (in the fourth row 
or column) have upper bounds ranging from 0.045 to 0.09. Tighter bounds on the geometric 
mean of two couplings are also found. In the Q = 2/3 case, the bounds are weaker. In that 
case, in fact, the \V t o\ and \Vjjb\ mixings are completely unconstrained, thus one could have 
large mixings between the third and fourth generations. The bounds on \ Vjjd\ an d \ Vu s \ are 
very weak, given by 0.15 and 0.56 respectively, while the bounds on \V u q\ and \V c q\ are as 
strong as the corresponding terms in the Q — —1/3 case. 

Other papers discuss certain particular processes in more detail, and look at constraints 
in more specific models. Let us first consider the case in which the isosinglet quark has 
Q = 2/3. This possibility was discussed in detail by Branco, Parada and Rebelo [ |156| . They 



showed that if one has an isosinglet Q = 2/3 quark, then the strongest signal will come from 
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D° — D° mixing. Suppose the Uu, Uc and Ut elements of the Q = 2/3 mass matrix are given 
by J u /Mu, Jc/Mu and Jt/Mjj respectively. They show that if one assumes that the J; are 
equal, then the current bound on D° — D° mixing gives an upper bound on J/Mu of 0.033. 
This bound could easily be saturated in realistic models, and thus D° — D° mixing gives the 
strongest constraint. If one were to take the Jj to be hierarchical, with a value of krrii, then 
the current limits give an upper bound on /c of 2, and thus the well-motivated k ~ 1 is well 
within reach of the next round of experiments. 

The Q = —1/3 case is more strongly motivated (since such states appear naturally 
in representations of Eg). These models have been analyzed in great detail in papers by 
Silverman and collaborators. The bounds on the Zds vertex coming from Kl — * A* + A<~> 
e, and the Ki — K$ mass difference were considered in early papers by Shin, Bander and 
Silverman [|159|1 and by Nir and Silverman ||160| , |161| ; bounds on the Zbd and Zbs vertices 
arising from B q — B q mixing and rare B decays were discussed in Refs. ||160|| and [|161|| , and 
followed up in Refs. ||162|| and [|163|| . Bhattacharya et al. ||164|| analyzed radiative B-decays 
in detail, and in another paper, Bhattacharya [|165|| analyzed the bounds from Z-decays. 
Although these processes are also discussed by BBP, they are described in much more detail 
in the above papers. More recently, Silverman [|166|| has analyzed the mixing constraints 
using the latest data from B physics, and looks at the constraints that will be reached in 
upcoming B factories. 

Lavoura and Silva [ I57 | extended the analysis to the case of multiple isosinglets. They 
also pointed out that a very strong bound comes from K + — > 7i + uV, a process which was 
not considered by BBP. As noted by Branco, Parada and Rebelo [ |156| |, using the realistic 
assumptions on the mixing mentioned above, the bound from K + — > 7T + uV is the most 
stringent for the Q = —1/3 case. The resulting bounds on the J/M are more stringent 
than in the Q = 2/3 case, J/M < .008. It is interesting that the rate for K + — > ir + vi7, 
which gives the strongest bound, has been measured and may be high (one event seen and 
a quarter of an event expected), but drawing conclusions on the basis of a single event is 
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certainly premature. 

The results of above two paragraphs apply to the case of an isosinglet vector-like fourth 
generation. The results are quite different in the case of a vector-like isodoublet. As shown 
in the analysis of the Aspon Model by Frampton and Ng ||168|| , the flavor changing Zdidj 
vertex will be suppressed relative to the isosinglet case by a factor of miirij /m 2 D . This is 
because, in the isodoublet case, the mismatch with light quarks occurs in the right handed 
sector, forcing a helicity flip of the usual quarks. (In the isosinglet case, the mismatch is 
in the left handed sector, thus no helicity flip is required.) This extra factor eliminates 
any significant constraints from flavor changing neutral currents. What about more exotic 
states? Recently, del Aguila, Aguilar-Saavedra and Miquel ||169|| looked at the constraints on 
anomalous top quark couplings in models with exotic quarks. They look at chiral and non- 
chiral singlets and doublets, including mirror quarks, and find some very general inequalities 
which allow one to go from LEP bounds on diagonal Z couplings to stringent bounds on the 
off-diagonal couplings. 

Thus, the bounds on mixing of a chiral fourth generation with the third are virtually non- 
existent, as are bounds on an isodoublet fourth generation, but the bounds on an isosinglet 
vector-like fourth generation are getting near the "interesting" range-and may be improved 
significantly with more measurements on K + — > ti + vV and on D° — D° mixing. 

It should also be pointed out that the four-generation model does have many additional 
phases. A detailed analysis of CP violation in the isosinglet Q = —1/3 case has been carried 
out by Silverman [|166j , |170|| . This is also discussed in BBP. The entire "unitarity quadrangle" 
is analyzed. We will discuss CP violation in more detail in Section VI. 

What are the theoretical expectations for the mixings? The Fritzsch ansatz (for the 



3x3 quark mass matrices) fails ||171| - |173|| for a 174 GeV top quark, although the generic 



expressions sin 8 = ^jm^/Mo or J nit/ My could easily be accommodated in other models. 

As noted above in the lepton case, many models with flavor symmetries will have the 2x2 

/0 A\ 

third-fourth generation mass sub-matrices of the form . These will have a 3-4 

\A B 
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mixing angle of 0(ym t /Mu), which is of order unity. Thus, one should keep in mind that 
the mixing angle between the third and fourth generation could be very large. As in the 
lepton case, one could imagine a symmetry in which the Q = 2/3 quarks are diagonal, and 
then the 3-4 mixing angle would be of order ^Jmb/Mo- 

The possibility that there is a symmetry prohibiting mixing altogether cannot be ex- 
cluded. With such a symmetry, the lighter of the U or D would be stable, leading to a 
cosmological disaster; however one could assume that Planck mass effects violate the sym- 
metry, giving a long (but possibly acceptable) lifetime of O(10 — 100) years. 

In the vector-like case, the mixing angles are also related to the J, discussed above, and 
the expectations are not very different. It should be noted that in the Aspon Model, the 
mixing angles are typically 10~ 5 — 10~ 3 in order to account for the appropriate amount of 
CP violation. 
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IV. LIFETIME AND DECAY MODES. 



In the previous Section, we discussed the masses and mixing angles of additional quarks 
and leptons. Now, we consider the lifetime and decay modes of such fermions. In the 
standard model, the lifetime and decay modes of the most recently discovered fermion, the 
top quark, were not particularly interesting-it was known that the top quark would decay 
very quickly (quickly enough that the width is large enough to obscure any structure in the 
toponium system) and that it would decay almost entirely into a b and a W. 

However, there are several interesting possibilities for the case of additional fermions. In 
the chiral case, the N could be heavier than the E, forcing the E to decay only via mixing; 
if the mixing angles are small (as discussed in the last Section), the lifetime could be very 
long. In the quark case, the mass of the D is likely less than the sum of the masses of the 
top quark and W, and thus the D will decay only via the doubly- Cabibbo suppressed c + W 
mode or the one-loop b+Z mode; either could give a long lifetime, especially if mixing angles 
are very small. In the non-chiral case, the GIM mechanism will not be operative, leading to 
tree-level flavor- changing neutral decays, such as E — > rZ, N — > u T Z and D — > bZ; these 
decays could give very unique and interesting phenomenological signatures. In addition, we 
will see that the mass-splitting of the leptons and quarks in the non-chiral doublet case is 
calculable, and gives lifetimes with potentially observable decay lengths. 

We will begin by discussing the lepton sector, first for the chiral case and then for the 
non-chiral case, and then turn to quarks. 
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A. Leptons 



1. Chiral Leptons 

Much of the early work on the phenomenology of heavy leptons ||174|| assumed that 
is less than mg, however, as discussed in the last Section, there is no particular reason for 
that assumption. Let us first assume the opposite-that m E < (and both greater than 
Mz/2. This case has been discussed in detail by Hou and Wong |[L75|| . The E will only be 
able to decay via mixing, E — > u T W*, where W* is a real or virtual W . The N will decay 
via either N — > EW* or iV — > tW*. The decay rates are 

T(N -> EW*) = 9 cos 2 6 34 (m 2 N /m 2 w , m 2 E /m 2 N )) , 

T(N -> tW*) = 9 sin 2 6 U (m^/m 2 ^, 0)) , (52) 

where f(a,/3) is given by ||176|| 

„, _ /(i-v^) 2 dx[(l - (3) 2 + (1 - g)x - 2x 2 ](l + (3 2 + x 2 - 2(f3 + f3x + x)) 1 / 2 

io [(l-x«)Hr«p] ( J 

This function accounts for both real and virtual W's. The rate for E — > v T W* is identical to 
the second of these equations with m N — > m E . Since the angle is expected to be small, one 
might expect that iV — > EW* would be favored over iV — > tW*, however, one must recall 
that the S and T bounds discussed above imply that the N and E must be fairly close in 
mass, and thus the decay might be significantly phase space suppressed. These rates are 
plotted in Figure 9. 
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FIG. 9. Relative branching ratios of the iV into EW* vs. tW* for various values of the mixing 
angle and the E to N mass ratio. One sees that the decay into tW* will dominate unless the 
mixing angle is very small, sin 2 26 is the mixing between the third and fourth generations. 

We see that the decay of iV into tW* tends to dominate, unless the mixing angle is 
extremely small. This leads to interesting phenomenology, as will be discussed in Section 
VII. If the E is heavier than the N, the results are the same with N «-» E and v T <-> r. 

Note that the decay rate of E — > u T W* is proportional to sin 2 6*34. In the previous Section, 
we noted that this angle could be very small. Simplifying the expression for the decay, and 
assuming that the mass of the E is greater than the W (thus the W is real), the width of 
the E is given by 

T(E -> v T W) = (180 MeV) sin 2 6 U ( — ) 3 (54) 

\m w J 

Consider the four plausible values of sin 2 6*34 discussed in the previous Section. If 
sin 2 #34 = m T /m E , then the decay is very rapid and would occur at the vertex. If 
sin 2 #34 = m^/mr?, then the lifetime is of the order of a few picoseconds, which has very 
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interesting phenomenological consequences — it might be possible to see the charged track. 
If sin 2 #34 = 0, then the E is stable; this would be cosmo logically disastrous, since they 
would bind with protons to form anomalously heavy hydrogen. If sin 2 #34 = m^/mp;, then 
the lifetime is approximately 10 — 100 years. We now address whether such a lifetime would 
pass cosmological muster. 

The bounds on the lifetime of a fourth generation charged lepton were first discussed 
several years ago ||177|| . They considered charged lepton masses ranging from 50 GeV to 



50 TeV. Stable leptons are ruled out by searches for heavy hydrogen. Since any decay of 
the E will result in photon emission, failure to observe such emission in the diffuse photon 
background implies that the lifetime must be less than 10 13 seconds (time of the cosmic 
background radiation (CMBR) production). Using COBE data on the CMBR, and requiring 
that the radiation in the decay not distort the CMBR, they found a bound on the lifetime 
which ranges from 10 9 — 10 11 seconds as the mass ranges up to 1 TeV. Very recently, an 
analysis by Holtmann, Kawasaki, Kohri and Moroi ||1 78|| looked at the radiative decay of 



a long-lived particle, X, and the effects on big-bang nucleosynthesis (the E — > u T W is not 
"radiative", however, a significant fraction of the energy will eventually turn into photons, 
and the results are the same). The photons emitted in the decay may photodissociate 
deuterium (lowering its abundance) and helium (which raises the deuterium abundance), 
destroying the agreement between theory and observation. They give bounds on the lifetime 
as a function of mxYx-, where Yx = nx/n^ is the relative abundance of the X. For heavy 



leptons, the abundance as a function of mass was calculated in Ref . ||177|| . For heavy lepton 
masses between 100 and 500 GeV, the contribution to flh 2 varies from 0.05 to 0.01, leading 
to a value of mxYx which varies from 6 x 10~ 10 to 1.2 x 10~ 10 GeV (for a Hubble constant 

one 



of 65 km/sec/Mpc). From Tables 3-5 of Holtmann, Kawasaki, Kohri and Moroi ||178 | 
can see that this correponds to an upper bound on the lifetime of between 10 7 and 10 8 
seconds. Given the uncertainties in the abundance calculation, nucleosynthesis calculation, 
deuterium and helium abundance observations, etc., this is not inconsistent with a lifetime 
of 10 — 100 years. Thus, the possibility that sin 2 #34 = m^/mp; is marginally allowed. 
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To summarize, if the mixing angle #34 is not very small, then both the N and E will 
decay via Cabibbo-suppressed decays: N — > tW* and E — > u T W*. If the angle is very small 
(of the order of m^ T /rriN or less), then the heavier of the two will decay into the lighter, 
while the lighter decays via the Cabibbo-suppressed decay. In this case, the latter lifetime 
could be quite long, as long as a few picoseconds for sin 2 9 U = m VT /m N and as long as 10 
years for sin 2 6*34 = mf^/m 2 ^. 



2. Nonchiral Leptons 

As noted in the previous Sections, an interesting feature of models with a vector-like 
doublet (with small mixing with light generations) is that the two members of the doublet 
are degenerate in mass, at tree level. This degeneracy will be split by radiative corrections, 
and the size of this splitting is crucial in understanding the lifetimes and decay modes of 
the heavy leptons. 



ZJ 



N 



FIG. 10. Diagrams contributing to the E-N mass difference. 



The splitting is due to the diagrams in Figure 10. The size of the splitting was first cal- 
culated by Dimopoulos, Tetradis, Esmailzadeh and Hall ||1 79|| (DTEH), and later calculated 
by Slier ||180|| (S) and even later by Thomas and Wells ||181|| (TW). The result is that the 
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charged lepton is heavier than the neutrino, with a mass splitting of 



a 



5m = -m z f(m 2 E /m 2 z ) 



(55) 



where 



/(x) = ^/ 1 C fo(2-z)lnfl+ X ). 
n Jo \ r(l — x) z J 



(56) 



For small x, f(x) — > 0, but for large x, f(x) — > 1, and thus the splitting reaches an asymptotic 
value of \amz — 350 MeV for m E » m Z - The splitting is plotted in Fig. 11. 
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FIG. 11. Mass difference between the E and the N as a function of the E mass. 



In DTEH ||179|| , the authors considered very heavy leptons (of the order of several TeV) 
and looked at the question of whether such leptons could constitute the dark matter (they 



can't). In S ||180|| , the decay width of E — > NeV and E — > N^V was calculated-the inverse 
of these widths corresponded to a lifetime of 1 — 2 nanoseconds, which is obviously of great 
phenomenological interest. In S, it was also pointed out that the value of the splitting is 
very robust, and that supersymmetric contributions to the splitting turn out to be very 
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small, since for most of parameter-space, the contributions turn out to be proportional to 
(\-sm 2 e w ). InTW plfl , it was pointed out that the dominant decay of the E will be into 
Ntt (a decay neglected in S), resulting in a much shorter (by roughly a factor of 10) lifetime. 
The decay length then is of the order of centimeters, rather than tens of centimeters, greatly 
complicating detection. In Fig. 12, we plot the decay distance in the lab frame as a function 
of rriE for a variety of center of mass energies. TW do propose an interesting signature 
involving triggering on an associated hard radiated photon; this will be discussed in Section 
VII. 

Thus, in the vector-like doublet case, without significant mixing with the lighter gener- 
ations, the charged member of the doublet is heavier and decays primarily into Ntt (with a 
VERY soft pion) with a decay length given in Fig. 12. 
20 




160 200 
M E (GeV) 

FIG. 12. Decay length for the E as a function of its mass. The lines are labelled with the center 
of mass energy of the collider in GeV. 



The fully leptonic decays, even though they have branching ratios of only a few percent, 
might be easier to detect, although even in that case, the very soft electron or muon would 
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be difficult to separate from backgrounds. If there is significant mixing, then both the E 
and the N will decay into the lighter generations. We now consider this possibility (which 
also applies to vector-like singlets). 

We now consider the case in which there is significant mixing of vector-like leptons. By 
"significant", we mean that the mixing angle, sin 2 # 34 , is greater than about 10~ n , so that 
the decay (see the expression for the lifetime above) occurs at the vertex. The results will 
be very similar to the chiral case, with one extremely important difference. Due to the 
breakdown of the GIM mechanism, the decays E — > tZ and N — > v T Z will occur. 

For rriE > Mz, the branching ratios of the E are given by ||182|| 



T{E -> tZ) 


U Et 


2 


(m| - 2m| + m|/m|)(m| - m 2 z ) 


T{E->v T W) 2cos 2 fl w | 


Uev t 


2 (m| - 2m 2 w + m%jm^){m 2 E - m^) 



(57) 



There is no particular reason to believe that this branching ratio should be small. One 
might expect \Uev t \ 2 to be of the order of m T /rriE, as discussed in detail in the last Section. 
As discussed in the last Section, in the isosinglet heavy lepton case, one finds that \Ue t \ is 
of order of the ratio of M 34 to M 44 in the leptonic mass matrix which would be similar to 
|Z7bi/ t |; t ne resulting branching ratio would be very large. In the isodoublet heavy lepton 
case, there is an additional suppression of m T /mE, which results in a branching ratio of 
about 0.1%. Thus, one expects a branching ratio for E — » tZ to be a fraction of a percent 
in the isodoublet case and very large in the isosinglet case. 

It is important to note that even if the branching ratio is as low as a fraction of a percent, 
the background for a particle decaying into tZ would be extremely small. A major problem 
with conventional heavy-lepton detection has been backgrounds; so the E — > tZ signature, 
even with a branching ratio of a fraction of a percent or so, might very well be easiest to 
detect. This will be discussed in more detail in Section VII. 
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B. Quarks 



The discussion of the lifetime and decay modes of the quarks follows the general discussion 
of those of leptons, with a few crucial differences. While the E and N are certainly much 
heavier than the r, it is not necessarily the case that the U and D quarks are much heavier 
than the top. In addition, the mass splitting in the vector-like quark case is a factor of three 
smaller than that for vector-like leptons, which can drastically affect the decay modes and 
lifetime. 



1. Chiral Quarks 

Due to constraints from the p parameter, the U and D quarks cannot have masses which 
are too different, and thus the decay U —>■ D + W or D —>■ U + W cannot occur into real 
Ws. Suppose, for the moment, that m\j > vtld. Then the allowed decays of the U will be 
U (D or g) + W*, U — > q + W, where W* refers to a virtual W. The allowed decays of 
the D are D — > (t or q) + (W* or W). In addition, one can have a flavor-changing neutral 
current decay D — > b + Z , which (in the chiral case) can occur through one loop. The fact 
that the flavor-changing neutral decay of a fourth generation quark could be significant was 
first pointed out by Barger, Phillips and Soni JTB3], and followed up by Hou and Stuart 
1. In the latter works, they noted that the decay D — ► b + Z dominates over other 



flavor- changing decays, such as D — ► b + 7 and D — > b + g. The possible D — > b + Z 
mode, like the E — ► r + Z mode, is very important phenomenologically due to the very clear 
signatures |[L86| , 187|| . Precise analytical formulae for the various decays, and a discussion of 



the decays of the D (from which much of the discussion below is taken), can be found in 
the recent work of Frampton and Hung [ |188|| . 

The two-body and three-body decay widths are given by the first of Eq. |52], with the 
obvious substitutions of |Vb g | or \Vjj q \ for cos 2 6 I 34 . The flavor-changing neutral current decay 
is given by 
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T(D -> bZ) = \V Dt 



G F {m D f 



{-^) 2 A(mu, m t ) I I 2 {m b /m D , m z /m D ) (58) 



4tt^2cos 2 0- 



14 ' 



where we have assumed, for simplicity, that Vub = —Vot so there will be a GIM suppres- 
sion when mjj = m t (this occurs in most models). We have also assumed that Vud is 
approximately unity Here, 



and rriheavy refers to the heavier of U and t. The factor I 2 is the standard two-body phase 
space factor, which is unity in the limit mo » rnz- 

The U decays very rapidly, at the vertex (unless the masses are unusually close together, 
which is unlikely in the chiral quark case, since they arise from different terms); however the 
D can be long-lived. The decay modes of the D depend crucially on the mass. 

First, suppose the D is lighter than the top quark. The decay can only occur into c + W 
or b + Z (one expects the decay into u + W to be highly suppressed). Two important issues 
arise: which of these has a larger branching ratio (since their phenomenological signatures 
are very different), and is it possible that the decay length might be large enough to be 
detected? We now address both of these questions. 




(59) 
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FIG. 13. Ratio of width of D -> b + Z to that of D -> c + W. Note the GIM suppression of 
the FCNC decay when the U mass equals the top mass. 

From the above formulae, the ratio of the widths can be calculated. The only dependence 
of the D-mass in this ratio is in the ratio of phase spaces-unless the D is fairly close in mass to 
the Z this ratio is near unity. The only mass-dependence is then in the U mass-dependence 
in A(mu, m t ) above. The result is given in Figure 13, for various values of the U mass (recall 
the fact that the U and D, in this range, cannot be different in mass by more than 20 GeV 
due to p parameter constraints). The results depend on the ratio of \Vd c \ to \ Vot\- Since the 
former is "doubly- Cabbibo-suppressed" (crossing two generations), one expects the ratio to 
be small. We see that the D — ► bZ decay mode dominates if the ratio is small, whereas 
D — > cW dominates if it is large (unless the U mass is very near the top mass. Since there 
is little theoretical guidance as to the size of this ratio, both signatures should be looked for. 

Could a vertex be seen? From the expression for T(D — > bZ), one can see that the width 



of this mode [ 184 1 is between 10 |Vbt| MeV and 10 |Vot| MeV over the mass range 
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of interest. This implies that the lifetime will be at most 6 x 10 _18 /|VDt| 2 seconds. This 
would mean that one must have |Vot| << 1CF 3 in order to detect a vertex. Although such 
a small angle is not expected in the chiral sequential quark case (the Aspon Model involves 
non-chiral quarks), it is not excluded, if one has a flavor symmetry with an almost unbroken 
3 + 1 structure. 

For rriD between 177 and 256 GeV, the D can decay into the top quark via the three-body 
decay into a virtual W. In this case, D — > t + W* will be competetive with D — > c + W 
and D — > b + Z — the latter are suppressed by the doubly-Cabibbo suppressed mixing angle 
or the extra loop, the former is suppressed by three body phase space. If we compare the 
rate for D — * t + W* with D — > b + Z, the mixing angle cancels out. In Figure 14, we have 
plotted the ratio of these two decays. 
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FIG. 14. Ratio of width of D -> t + W* to that of D -> b + Z for various values of M v . The 
non-chiral line corresponds to the vectorlike doublet case, which is independent of Mjj. For the 
non-chiral isosinglet case, the ratio is too small to be seen on the graph. 
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Thus, as pointed out in Ref. ||188|| , the three-body decay of the D is irrelevant for D 
masses below about 230 GeV (depending on Mu), and thus the arguments of the previous 
two paragraphs are unchanged. 

As the mass increases beyond 230 GeV, the three body decay becomes more important, 
and as soon as the mass exceeds 256 GeV, the two-body decay D — > t+W becomes accessible. 
At that point, D —* t + W becomes the dominant decay. Again, very small mixing would 
be needed to see the vertex, \Vm\ < 10~ 5 . 

What about the decay of the U quark (again, assuming that mu > iud)? The U can 
decay into D and a virtual W, U — > D + W*, or into a light quark (most likely a b) and 
a real W, U — > b + W. Which decay mode dominates depends on how close the U and D 



are in mass and on the \Vub\ mixing angle. This was discussed in Ref. ||188|| . As nijj varies 
from 180 to 250 GeV, the width for U -> b + W varies from l.7h\V Ub \ 2 GeV to A.l\V Ub \ 2 
GeV. For U — > D + W*, the width (assuming \Vud\ is nearly unity) is 5.2 x 10~ 5 GeV for 
n^u/ m D — 1-1, and drops to 3 x 10~ 8 GeV for mu/m^, = 1.02. For moderate mixing angles, 
\Vjjb\ > 0.003, the [/—>& + W decay mode will always dominate. For mixing angles in 
the range between 10~ 3 and 10~ 4 , which dominates dependes sensitively on the mass ratio 
between the U and the D. For mixing angles below 10~ 4 , the U — ► D + W* will dominate. 
In no cases will the width be small enough so that a vertex can be seen: the U will decay 
at the vertex. 

If run > mu, everything we have said above will carry through, exchanging t with b, etc. 
A principal difference is that one can consider lighter long-lived U quarks (since one can 
have mu < m t ) than for D quarks. 

To summarize, if mu > mo, then the primary decays modes of the D will be D — >• c + W 
and D — > b+Z if the D mass is below about 230 GeV. The former will dominate if | V^d|/| Vut\ 
is greater than 0.01, the latter will dominate if it is less than 0.001. As the D mass increases 
above 230 GeV, the decay D — ► t + W* begins to dominate. In all cases, one must have very 
small mixing angles (\Vot\ < 10~ 3 or 10~ 5 , depending on the mode) in order for the decay 
to occur a measureable distance from the vertex. The U will always decay at the vertex, 
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into either D + W* or b + W, depending on the precise masses and mixing angles. 



2. Nonchiral Quarks 

For vector-like quarks, either isosinglet or isodoublet, many of the above results are 



unchanged. As noted in Ref. ||188|| , the decay widths into real or virtual W's will not be 
significantly changed. As in the nonchiral lepton case, there are two major differences: the 
mass difference between the U and the D in the isodoublet case is calculable (if mixing is 
small), and the flavor- changing neutral decay D — > bZ can occur at tree level. 



In the isodoublet case, the mass difference between the U and the D can be shown ||180 
to be 1/3 that of the lepton case, and will thus be between 70 and 110 MeV, if the mixing 
with lighter generations is small. This means that the hadronic decay is forbidden, and the 
only decay would be U — > Dei/, with a lifetime of the order of microseconds. Thus, in the 
absence of mixing, both the U and the D would be absolutely stable, as far as accelerators 
are concerned, and the mass difference is irrelevant. We thus must consider both U and D 
decays. 

The D can decay, as in the chiral case, into either c + W, t + W*, t + W or b + Z. The 
only difference in the above discussion is the rate for D — > b + Z. In the expression for 
T(D — > b + Z) in the last section, one must replace the factor of g 4 A(mu, m t ) / 647r 4 from the 
loop with ml/m 2 D , in the isodoublet case, and with unity in the isosinglet case (the factor 
of \VDt\ 2 is the same in each case). In the isodoublet case, this does not change the result 
for the width by more than an order of magnitude, but in the isosinglet case, it increases 
the width by roughly three orders of magnitude. 

First, if the D is lighter than the top quark, then the decay modes in the isodoublet 
case are very similar to the chiral case — the results depend sensitively on the ratio of \Vd c 
to |Vbt|, and either the D —>■ c + W and D — > b + Z will be dominant. A displaced vertex 
can only be seen if |Vbt| < 10~ 3 . This is expected in the Aspon Model case, and thus one 
can expect a significantly displaced vertex in the model (as emphasized in Ref. |]188|| ). In 
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the isosinglet case, the absence of the m1/m 2 D suppression makes the D — > b + Z mode 
dominant, and also requires that \Vr>t\ « 10~ 5 in order for a displaced vertex to be seen 
(which is not expected in the Aspon Model). 

If the D is heavier than the top quark, then the argument in the previous section carries 
through without significant modification in the isodoublet case; the cross-over where the 
three body decay becomes relevant is closer to 210 GeV than to 230 GeV. In the isosinglet 
case, however, the D — > b + Z mode dominates, even over the two-body top quark decay, 
until m_o is well over 300 GeV, and remains significant even to much higher masses. 

Thus, we see that in either the chiral or non-chiral case, the decay D — > b + Z is very 



important. It was pointed out in Ref. [|188|| that, if this mode were detected, then one could 
look at the chirality of the Z to determine which of the two cases applies. This might be 
the quickest way to determine the chirality of the heavy quarks. 

Unless the U is very heavy, above 310 GeV, the flavor- changing neutral decay U — > t + Z 
is forbidden or highly suppressed by phase space. Since U — > c + Z is suppressed by small 
mixing angles, one has the U — > b + W decay mode dominating. Again, one can only detect 
the vertex if the angle is very small, less than 10~ 5 . 
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V. DYNAMICAL SYMMETRY BREAKING 



It is fair to say that perhaps one of the most important discoveries that can be made 
in the future will be that of the Higgs boson. In the absence of any alternative plausible 
explanation for particle masses, the concept of spontaneous breakdown of a gauge symmetry 
via the Higgs field as the origin of all masses is universally accepted. The only problem is 
that it has not been found. Not only does one not know its mass, but one also does not know 
in what shape or form it should be. One fact that we do know however , regardless of how 
massive and in what form the Higgs boson may be, is the scale of electroweak symmetry 
breaking: v = 246 GeV. All masses in the SM are expressed in terms of that scale, e.g. 
Mw = (1/2) gv and m« = gY^v/y/2, where gy i are Yukawa couplings. There is a rather 
intriguing fact: With v/y/2 ~ 174 GeV, it follows that the top quark Yukawa coupling, g t , 
is of order unity, unlike all other fermions. That the top quark is so heavy and its mass is 
so close to the electroweak breaking scale is cause to wonder about any relationship that it 
might have with the mechanism of symmetry breaking itself. 

If the top quark mass is so close to the electroweak scale, is it possible that it itself is re- 
ponsible for the breaking of the electroweak symmetry? This fascinating possibility was first 
entertained by Refs. P5| - [f46f . It is now commonly refered to as "top-condensate models". 
In these models, the Higgs boson is generally viewed as a ti composite field generated by 
some unknown dynamics at a scale A ^> 246 GeV. In other words, the electroweak symmetry 
breaking is dynamical in these scenarios, i.e. it is broken by a ti condensate or something 
similar. One interesting feature of these models is the prediction of a heavy fermion mass 
(e.g. the top quark mass or a fourth generation mass) as a function of the Higgs boson mass 
through the so-called compositeness conditions. Before describing these models and their 
variants, let us first summarize the results of the simplest version, that of Bardeen, Hill and 



Lindner |[46|| . In Ref. p6| , the first scenario including only a heavy top quark predicted a top 
mass of order 230 GeV. This is now excluded. Ref. |46| also presented results involving a 
heavy, degenerate fourth generation quark doublet, assuming that the top quark and other 
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fermions are much lighter. This result combined with the known top quark mass indicates 
that something else other than the top quark alone must exist if this picture has any chance 
of being correct. 

The subject of dynamical symmetry breaking a la top-condensate deserves a whole re- 
view; a book on the topic already exists |189 |. The best we could do here is to review salient 



points and results, especially those pertaining to the subject of this review. We first state 
the so-called compositeness conditions used in the top-condensate type of models. As we 
shall see below, these conditions allow us to relate the masses of the heavy fermions to that 
of the Higgs boson. Simply speaking, it is the requirement that the Higgs quartic coupling 
A and the relevant Yukawa couplings diverge at the same scale Ac (the Landau poles) while 
the ratio of A to the square of the Yukawa coupling(s) remains finite, namely 

A(/i), g t {n) ^ oo (60) 

*(aO/s?(aO m —^ c const - ( 61 ) 

where gt can be either the top or just a generic Yukawa coupling. The above boundary 
conditions modify the structure of the SM Lagrangian at the scale Ac in the following 
way. Let us consider a simple toy model where there is a degenerate heavy quark doublet 
Q = (U,D) coupled to the Higgs field $ (light fermions will be ignored in this discussion). 
Let us rescale $ as follows 

$ — > %/g f , (62) 

where gj is the Yukawa coupling of the degenerate quark doublet. The SM Lagrangian (with 
only that degenerate doublet present) becomes 

C = C kinetic (U, D) + ZtD^lD^o + m 2 $J$ - ^(<4$ ) 2 + Ql%D r + Q L $%U R + h.c, 



(63) 



where 



$^ = ia 2 %, = l/gj, m 2 = Z^m 2 , and A = Z%\. (64) 
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The first remark one can make when one looks at the above expressions is the composite- 
ness condition itself: the vanishing of the wave function renormalization constant. Indeed, 
as can be seen from the expression for one immediately notices that, if gf has a Landau 
singularity at some scale Ac, the wave function renormalization constant for the Higgs field, 
Z\phi, vanishes at that same scale. In other words, the Landau pole is identified with the 
compositeness scale. Furthermore, if A and gf develop a singularity at the same scale in 



such a way that conditions (|60[) are satisfied, the Lagrangian at the compositeness scale 
becomes 

C = CkinetUU, D) + Q L %D R + Q L ^U R + m 2 $ $ + h.c. (65) 

$o is now just an auxiliary field and can be integrated out, resulting in a Nambu-Jona-Lasinio 
form for the Lagrangian, namely 

C = C kinetic (U, D) + G Q L (U R U R + D R D R )Q L , (66) 

where Go — —1/m 2 . In this picture, the Higgs boson becomes a fermion-antifermion com- 
posite particle below the scale Ac- 

What might be even more interesting is the relationship between the Higgs boson mass 
and that of the heavy fermions in the top-condensate type of scenario. In particular, the 
Higgs mass, m#, can be seen to be bounded from above by 2m j and from below by rrif, 
where mj is a heavy fermion mass. The search for a heavy fermion is intimately tied to 
the search for the Higgs boson-a feature already seen in the discussion of gauge coupling 



unification [65]. To see this in the context of the top-condensate type of model, let us look 
at the RG equations for A and gf, the heavy fermion Yukawa coupling, at one loop level 
(the one- loop terms of Eq. ( |47a| )) in the toy model with one doublet of heavy quarks. Let 
us define x = A/g 2 . The one-loop RG for x is then 

16vr 2 ^ = 24#J(x-x + )(x-x_), (67) 

where x± = (— 1 ± 3)/4, if both members of the quark doublet are degenerate in mass, 
or i± = jq( — 1 ± V65), if one member is much heavier than the other one (e.g. the 3rd 
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generation case). The boundary conditions fl60"D are satisfied if x is one of the two fixed 
points, x±. The solution x_ (always negative) is ruled out by vacuum stability. This leaves 
us with x + . From the definition of mj and m#, one can then obtain a relationship between 
these two masses as follows: 



m 2 H = Am 2 f x + . (68) 



In the large N c limit, the right hand side of Eq. (|67| ) becomes 4N c (x — l) with the fixed point 
being x = 1 which implies m# = 2m f, a familiar result found in the Nambu-Jona-Lasinio 
model. Since x + is always greater than 1/4, one can easily see that, for finite N c , one obtains 
the bound 

rrif < m H < 2rrif. (69) 

As we have mentioned above, a necessary condition for this scenario to work is the 
boundary (J§Q). The minimal SM with three generations unfortunately does not satisfy 
these boundary conditions: the top quark with a mass of 175 GeV is simply too "light". 
One needs eiher a heavier quark (216-230 GeV), which cannot be the case with the minimal 
SM, or more heavy quarks or leptons such as in the four generation scenario. The four 
generation case was studied by Hill, Luty and Paschos |190| in the context of Majorana 



neutrinos. Later, Hung and Isidori [|117|| also examined the four generation scenario as 
part of an overall "anatomy" of the Higgs mass spectrum, starting form m# of 65 GeV to 
mjj > 2mz- It was in this last mass range, m# > 2mz, that the focus on top-condensate 
types of models was diverted to. In their analysis, Ref. ||117|| assume a Dirac mass for the 



fourth neutrino and degenerate quarks and degenerate leptons for the fourth generation. It 
turns out that the inclusion of such a fourth generation drastically modifies the evolution of 
the couplings even if these fermions were lighter than the top quark. In the analysis of Ref. 
117|| , Eqs. ( |47a| ) were used at the one loop level. The results are shown in Fig. 15 below, 



where the mass of the top quark is fixed at 175 GeV and that of the fourth lepton doublet 
is fixed at 90 GeV. 
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FIG. 15. RG evolution of the Yukawa couplings gf (full and dash-dotted lines) and g 2 q (dashed 
lines). The top mass is always fixed to be 175 GeV and the heavy-lepton mass is assumed to be 
90 GeV. The dash-dotted line is the evolution of g 2 without extra fermions. Near each dashed 
line is indicated the value of m q and the corresponding value of uih obtained by the requirement 
Al = A q (the error on both mg and m q is about ±10 GeV. 

A few remarks are in order concerning the above figure. As the figure caption already 
indicates, the dash-dotted line represents the top Yukawa coupling squared, g 2 , as a function 
of energy for the minimal SM. One can see that g 2 actually decreases in value with energy 
and remains finite at the Planck scale. This is so because a mass of 175 GeV is "too small" 
to provide a large enough initial value plus the fact that the contribution from the QCD 
coupling to the (3 function for g 2 occurs with the opposite sign. When a fourth generation is 
added, the evolution of the couplings change drastically, partly due to the fact that there are 
more degrees of freedom than the minimal case. So, as long as the initial Yukawa couplings 
of the fourth generation are not too small, there will be Landau poles that appear below 
the Planck scale. As one can see from Fig. 15, this occurs when the fourth generation quark 
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mass ttiq > 150 GeV. Also an interesting phenomenon occurs: all Yukawa couplings "drag" 
each other in such a way that they all "blow up" at the same point. The top quark Yukawa 
coupling, which by itself decreases with energy, is now "dragged" in such a way as to develop 
a Landau pole at the same time as the fourth generation quarks. These results- the existence 
of Landau poles below the Planck scale in the presence of a fourth generation- encourage a 
reconsideration of the top-condensate type of model. 

Since the quartic coupling A is a free parameter, one can now choose its initial value, 
i.e. its mass, in such a way that the boundary conditions (|60|) are satisfied. This procedure 
results in a relationship between the fourth generation quark mass and the Higgs mass. This 
is shown in Fig. 15 (m,f is the notation used in the figure to refer to the fourth generation 
quark mass instead of mq which is used here) . The relationship for low 4th generation quark 
mass (e.g. 150 GeV) is strikingly similar to the one obtained in the discussion of the gauge 



coupling unification of Ref. [p5[ . The range of mass used in Ref. ||117|| for ttlq is however 
considerably larger than that of Ref. |H| because, in that study of a top-condensate type of 
model, no constraint on where the Landau poles should be has been used. In fact, for masses 
larger than 150 GeV, the Landau poles move down in energy. For example, if mq = 230 
GeV which corresponds to m# = 300 GeV, the Landau pole is situated at approximately 



100 TeV. As remarked by Ref. 1|117| , the relationship between the Higgs mass and the fourth 



generation quark mass approaches more and more the fixed point value for a degenerate 
quark doublet, ran — V^mQ, as the Landau pole "approaches" the electroweak scale. This 
can be seen in the figure shown above. 

The scenario described above intimately links the search for the fourth generation quarks 
and leptons to that of the Higgs boson, and vice versa. It is not clear at this point how the 
higher mass values (> 160 GeV) for the fourth generation quarks would affect the evolution 
of the gauge couplings, except for the fact that perturbation theory breaks down above the 
Landau poles and it is not legitimate to evolve those couplings beyond that point. For 
the purpose of this report, we shall however leave open the possibility of such top-quark 
condensate type of models as a possible mechanism for electroweak breaking. It is partly 
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for this reason that the range of masses considered in the search for long lived quarks by 



Frampton and Hung ||188j| is larger than the mass range used in the study of gauge coupling 



unification of Hung [55 



Recently |[L91]| , there was an analysis of the two-loop RG equations in the SM with three 
and four generations. Although many results presented there ||191|| were already discussed 



in [65| and ||1 1 7|1 , there were statements that are not correct. In order to clarify the issues, 



we repeat the statements of Ref. |191| and explain why they are misleading. Ref. |191| chose 



the masses of the leptons and the quarks of the fourth generation as follows: m^jniQ = 1/2 
and 1, and restricted tuq to be greater than 180 GeV which was referred to by the authors 
as the direct experimental constraint. Furthermore, an upper bound of 200 GeV for nin 



was used. Using these constraints, Ref. [|191|| claimed that a fourth generation is ruled out 



by plotting the allowed regions in the m# — mq plane (Fig. 11 of ||191|| ). First, tuq less than 
180 GeV is not ruled out by direct experiment if a long-lived quark decays in the detector 



at a distance between 100 /zm and 1 m, a subject discussed by Frampton and Hung ||188 
As we shall discuss in the section on experiments, there is and will be such a search at the 
Tevatron. Second, the constraints on at the present time from precision experiments 



are rather loose at best. A much larger bound than 200 GeV is possible |192| . For example 



Refs. |193| - |195|1 gave an upper bound as high as ~ 280 GeV, while Ref. [p.96|l gave an upper 



bound (within 95 % CL) from 340 GeV to 1 TeV. In summary, the mass ranges used by 



191|| to rule out a fourth generation are not warranted. In fact, the values used by Hung 
|65| , itiq = 151 GeV and m# = 188 GeV, were shown to lead to a better unification of the 
SM gauge couplings, and higher masses (such as 180 GeV) for rriq were not used because 
precisely of the fact that the Landau poles were much too low to trust perturbation theory 
in the evolution of the gauge couplings. 

As the above discussion and earlier ones made it crystal clear, there are several theoretical 
reasons for looking at quarks and leptons beyond the third generation. In particular, our 
primary motivation in this review is to examine those reasons which "predict" fermion 
masses which are within reach either of present experiments such as the ones performed 
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at the Tevatron or at future machines such as the LHC, NLC, etc. By "within reach", we 
mean masses which are close to the top quark mass ( ~ 150 -230 GeV) for the quarks. It 
goes without saying that there exists also plenty of reasons for considering fermions which 
are much heavier than the top quark, some of which are even primary focus of the topics 
discussed in this section: the techniquarks and leptons of Technicolor models [|1971-|200 



It is outside of our goal to review such topics but for completeness, we shall give a brief 
description of one of such scenarios: the topcolor assisted Technicolor model. 

Our discussion of top-condensate type of models above relies on one crucial assumption: 
The entire electroweak symmetry breaking is due to the condensate of the top quark and/or 
that of the fourth generation quark. The topcolor assisted Technicolor model of Hill ||201|| 



did away with that assumption. But then how does one explain the fact that the top 
quark mass is so much larger than all other fermion masses and so close to the electroweak 
breaking scale? The salient points of the topcolor model are basically the assumptions that 
the electroweak symmetry is still broken by some form of Technicolor, there exists an extra 
("topcolor") group, S77(3)i ® U{1)y 11 which couples preferentially to the third family and 
which triggers a dynamical condensate for the top quark, giving rise to a large top mass. 
This new "topcolor" group is assumed to be spontaneously broken by Technicolor at a scale 
~ 1 TeV. An extended Technicolor interaction (ETC) is also assumed so that quarks and 
leptons can obtain some mass which is much smaller than the top mass. The top quark 
itself obtains most of its mass from the condensate with a small contribution coming from 
ETC. Variations of the model include cases in which there are SU(2) singlet quarks which 
are however very massive (~ 1 TeV). In any case, all fermions beyond the 3rd generation is 
these models have masses around 1 TeV. Although there is a rich phenomenology involving 
objects such as top-pions, etc., it is beyond the scope of this review to discuss it here. 

We end this section by mentioning that there are interesting recent developments on the 
subject of electroweak symmetry breaking involving the so-called seesaw mechanism of quark 
condensation ||202| , |203|| . This is yet another variant of the topcolor model. In this particular 
scheme, the top quark mass obtains its "observed" value dynamically by a mass mixing with 
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a new SU(2) singlet quark x °f the form fi x tXLtn + fnt x tLXR + HxxXlXr with m <x ~ 0.6 
TeV, /x x< ~ 0.9 TeV and fi xx ~ 2 TeV. The physical top mass is found by diagonalizing that 
2x2 matrix giving m t as ji xt m tx / /i xx . The other mass eigenstate would be ~ 2 TeV and it 
might not be easy to be observed directly. The model predicts a number of pseudo Nambu- 
Goldstone bosons, some of which are x~ bound states. How this new degree of freedom 
manifest itself experimentally is a subject which is under active investigation. 

One last comment which is worth mentioning here is the mass of the Higgs boson in the 
different variants of the topcolor model. Generically, the physical Higgs scalar would have 
a mass of the order of a TeV just like standard Technicolor models. However, in the seesaw 
version of the topcolor model, there appears an extra "light" Higgs scalar (in addition to 
the charged ones) as one would have in a two-Higgs doublet model. Depending on a delicate 
cancellation of some parameters of the model, one could have a "truly light" scalar with 
mass of 0(100) GeV. These remarks are meant to emphasize the importance of the search 
of Higgs scalar(s) in addition to that of new quarks and leptons. 



VI. CP VIOLATION. 
A. CP Violation in the Standard Model. 

At first sight, there may appear to be no connection between additional quarks and lep- 
tons and the violation of CP symmetry. The object of this subsection is therefore to convince 
the reader that the better understanding of CP violation may necessitate the incorporation 
of additional fermions. This is not an inevitable conclusion but it is a suggestive one. 

To set the scene, we need to describe the status of CP symmetry and its violation in the 
context of the unadorned standard model. This discussion will be in two separate parts: 
weak CP violation (the KM mechanism), and the strong CP problem. The former is not 
really a problem for the standard model, merely one that is not yet verified unambiguously 
by experiment. The latter, the strong CP problem, is definitely a difficulty, a shortcoming, 
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for the standard model and one whose solution (we shall mention the axion possibility) is 
still unknown. 

The gauge group of the standard model is SU(3)c x SU {2)l x U(1)y , broken at the weak 
scale to SU(3)c x U(l)y- Under the standard group the first generation transforms as: 



Qi 



yd j 



5 °L 



,u L ,d L ; L L 

and the second (c, s, z/ M , /x) and third (i, b, u T , r) generations are assigned similarly. 



(70) 



The quarks acquire mass from the vacuum expectation value (VEV) of a complex SU (2) L 



doublet of scalars 



giving rise to up and down quark mass matrices: 



M(U) = Ag < cfP >■ M(D) = \° < cfP > 



(71) 



which are arbitrary matrices that may, without loss of generality, be chosen to be hermitian. 
The matrices M(U),M(D) of Eq. (|7l|) are defined so that the Yukawa terms give e.g. 
QlM(U)ur + h.c. and can be diagonalized by a bi-unitary transformation: 



K{U) L M(U)K{U)l l =dtag(m U} m C} m t ) 



(72) 



K(D) L M(D)K(D) R 1 = dtag(m d ,m S} m b ) 



(73) 



These mass eigenstates do not coincide with the gauge eigenstates of Eq. (|70|) and hence the 
charged W couple to the left-handed mass eigenstates through the 3x3 CKM matrix Vqkm 
defined by: 



Vckm = K{U)lK{D)^ j 1 



(74) 



This is a 3 x 3 unitary matrix which would in general have 9 real parameters. However, 
the five relative phases of the 6 quark flavors can be removed to leave just 4 parameters 
comprising 3 mixing angles and a phase. This KM phase underlies the KM mechanism of 
CP violation. 



With N generations and hence a N x N mixing matrix there are N(N — l)/2 mixing 
angles and (N — l) 2 parameters in the generalized CKM matrix. The number of CP violating 
phases is therefore (N - l) 2 - N(N - l)/2 = (JV - 1)(N - 2)/2. This is zero for N = 2, 



one for N = 3, three for iV = 4, and so on. In particular, as Kobayashi and Maskawa |27 
pointed out, with three generations there is automatically this source of CP violation arising 
from the 3x3 mixing matrix. This is the most conservative approach to CP violation. This 
source of CP violation is necessarily present in the standard model; the only question is 
whether it is the only source of CP violation. Since the only observation of CP violation 
remains in the neutral kaon system, there is not yet sufficient experimental data to answer 
this question definitively. 

There are various equivalent ways of parametrizing the CKM matrix. That proposed 
27] by KM involved writing: 

/ . „ „ . „ . „ \ 



CKM 



cosOi —sinOicosOs —sin9\sin9 3 

sin0icos92 cos9iCos92Cos9 3 — sin92sin9 3 e lS cos9\Cos92sin9 3 + sin92Cos9^e lS 
sin9\sin92 cos9isin9 2 cos9 3 + cos92sin9 3 e tS cos9isin92sin9 3 — cos9 2 cos9 3 e lS 



(75) 



Another useful parametrization |[204| writes: 



V, 



CKM 



|A 2 



(76) 



A X 3 A(p - it]) 
-A 1 - |A 2 X 2 A 

\ 3 A(l-p-r]) -X 2 A 1 

In Eq.(TS), A is the sine of the Cabibbo angle sind\ in Eq.([75]) and CP violation is propor- 
tional to 7]. If we write the CKM matrix a third time as: 



V 



V u d V us V u b ^ 



V, 



CKM 



then the unitarity equation (VckmYVc 



V c d Vcs V c b 

Vtd V ts Vtb 
ckm = 1 dictates, for example, that 



(77) 



/ 
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V* b V ud + V; b V cd + V: b V td = (78) 

This relation is conveniently represented as the addition of three Argand vectors to zero in a 
unitarity triangle. Dividing out the middle term of Eq. (|78| ) and using the parametrization of 
Eq. ([75]) leads to the prediction of the standard model with KM mechanism that the vertices 



of the unitarity triangle in a p — 77 plot are at the origin (0, 0), at (1, 0) and (p, rf). Thus, the 
area of the unitarity triangle is proportional to 77 and hence to the amount of CP violation. 

The measurement of the angles and sides of this unitarity triangle are the principal goals 
of the B Factories (see e.g. ||205|| for a review)^. As we shall see, alternative models will 



give quite different predictions and hence be easily distinguishable from the KM mechanism 
when accurate measurements on B-meson decays are made in B Factories. 



B. Strong CP and the Standard Model. 

Next we turn to a brief outline of the strong CP problem in the standard model. (More 
detailed reviews are available in ||207| - [209|j ). The starting observation is that one may add 



to the QCD lagrangian an extra term: 

L = Y1 Qk(hn D n ~ m )<?fc - ®G^G^ (79) 

k 

where the sum over k is for the quark flavors and is the partial derivative for gauged 
color SU(3)c- The additional term proportional to violates P and CP symmetries. This 
term is a total divergence of a gauge non-invariant current but can contribute because of 
the existence of classical instanton solutions. It turns out that chiral transformations can 
change the value of G via the color anomaly but cannot change the combination: 

6 = 6- argdetM(U) - argdetM(D) (80) 



A recent prelimary report from CDF [20£] gives a value of 0.8 ± 0.4 for sin 2(5. 
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where det M(U, D) are the determinants of the up, down quark mass matrices respectively. 
Thus which is an invariant under chiral transformations measures the violation of CP 
symmetry by strong interactions. A severe constraint on arises from the neutron electric 
dipole moment d n which has been measured to obey d n < 10~ 25 e. — cm. ||210| , [211|| . A 
calculation of d n | |212| , |213| | leads to an estimate that < 1CT 10 . This fine-tuning of is 
unexplained by the unadorned standard model and raises a serious difficulty thereto. 

A popular approach (which does not necessitate additional fermions) involves the axion 
mechanism which we briefly describe, although since only a relatively narrow window remains 
for the axion mass, and since the mechanism is non-unique, it is well worth looking for 
alternatives to the axion for solving the strong CP problem. 

In the axion approach, one introduces a color-anomalous global U(l) Peccei-Quinn sym- 
metry ||214| , pi5|| such that different Higgs doublets couple to the up- and down- type quarks. 



The effective potenetial now becomes a function of the two Higgs fields and 0(x) regarded 
as a dynamical variable. An analysis then shows that the potential acquires the form: 

V = V(H h H 2 ) - cosG (81) 

and hence the minimum energy condition relaxes to zero. 

Because a continuous global symmetry is spontaneously broken, there is a pseudo- 



Goldstone boson ||216| , |217|1 , the axion, which acquires a mass through the color anomaly 
and instanton effects. The simplest model predicts an axion with mass of a few times 
lOOkeV, but this particle was ruled out phenomenologically. Extensions of the axion model 



[ 218H221| lead to an axion mass which becomes a free parameter. Empirics constrain the 



mass to lie between about a micro-electron Volt and a milli-electronVolt, and searches are 
underway for such an axion. 



In one variant of the axion approach, new heavy quarks are necessary |[218[ |; in fact, this 
was the first-proposed model to avoid the experimentally-excluded "visible" axion in favor 
of an "invisible" axion. Alternative versions of the invisible axion [ |220| , |221|| do not involve 
extra quarks. 
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A second solution of the strong CP problem is to assume that m u = 0, but this seems 
to be at variance with the successes of chiral perturbation theory ||222 - 225| . 



C. Strong CP and Extra Quarks. 

For present purposes, let us believe neither the axion nor the massless up quark. In this 
case we are inevitably led to the existence of further fermions beyond the standard model. 
Such a conclusion can have the additional bonus of connecting the strong CP solution to 
the occurrence of weak CP violation. 

The appearance of new fermions in this context was first suggested by Nelson [ |22(j| | and 
Barr [227 , 22j| who looked within the framework of GUTs. Their additional fermions have 



gigantic masses ~ 10 12 GeV, well beyond accessible energy, but their basic idea involving the 
quark mass matrix texture is one that will reappear in the models where the new fermions 
are at accessible masses. 



Nelson |[226| invented a GUT based on SU(5) gauge x (SO(3) x U(l)) g i b a i- The fermions 
are assigned to the representations: 

[10.3.0] + [5, 3, 0] + [10, 1, 1] + [10, 1, -1] + [5, 1, 1] + [5, 1, -1] (82) 

Note, in particular, the additional fermion representations in the last four terms of Eq. flS2]) . 

The scalars are the (5, 1,0) which is complex and contains the standard Higgs doublet 
together with a superheavy color triplet; then there are (rj,3,0) where are SU(5) repre- 
sentations containing singlets of SU(3) x SU{2) x U{1) such as 1, 24, 75. 

If we write the most general Yukawa terms we find that the couplings of the quarks to 
the light Higgs is complex but has a real determinant. Thus, if the lagrangian respects CP 
symmetry, the value of vanishes at tree level. 

Looking at the loop corrections to G, it turns out that the additional fermions must be 
lighter than the GUT scale by three or four orders of magnitude to suppress adequately. 

The secret of the Nelson model lies in the arrangement of additional fermions such that a 
basis may be chosen where complex entries in the quark mass matrix are always multiplied 
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by zero in the evaluation of the determinant. Barr |[227| , [228]1 examined what are the general 



circumstances under which this suppression of happens. 

Barr was led to the following rules for a GUT in which = Qqcd + ®qfd = at 
tree level. Let the gauge group be G and let CP be a symmetry of the lagrangian. Then 
0QCD — and the couplings have no CP violating phases. Let the fermion representation 
be divided into two sets F and R where F contains the fermions of the three families and R 
is a real non-chiral set. Let R be composed of a set C and its conjugate set C. Then the 
following two conditions are sufficient to ensure that = at tree level: 

• At tree level there are no Yukawa mass terms coupling F fermions to C fermions, or 
C fermions to C fermions. 

• The CP violating phases appear at tree level only in those Yukawa terms which couple 
F fermions to R = C + C fermions. 

In such Nelson-Barr GUT models, strong CP is solved by physics (the additional 
fermions) close to the GUT scale. The unrelated (in this approach) physics of weak CP 
violation arises from the usual KM mechanism. 

Non-GUT models which adopt the Nelson-Barr mechanism have been discussed especially 
in the papers ||229| p33f| . The model proposed in ||229|| by Bento, Branco and Parada intro- 



duces a non-chiral charge —1/3 quark together with a complex singlet scalar S. The field S 
develops a complex VEV < S >= Ve ia while the standard Higgs doublet has VEV < >= v 
which is real. The KM phase 8km is generated from a in an unsuppressed manner. is zero 
at tree level and its loop corrections are suppressed by powers of (< > / < S >) = {v /V). 
Consequently can be naturally sufficiently small; e.g. if V > lOOTeV, the Yukawa can be 
~ 1CT 1 while even if V > ITeV, the Yukawa coupling can still be as large as ~ 10~ 2 . Further 



papers examine other consequences of such a model. In ||230|| the impact foiB — B mixing 
is found, while in [231] D — D is also found to differ from the standard model predictions. 



One general feature of this type of model is that is suppressed by the Nelson-Barr type 



86 



mechanism and, as in Nelson-Barr, the CP violation arises from 8 km- In the Aspon Model 
described below, the CP violation necessarily arises from an additional mechanism because 
there the generation of 5km is highly suppressed. 



D. CP and a Fourth Generation. 

Before moving to that, let us mention a number of valuable papers which discuss the 
parametrisation of the CKM matrix when the number of generations of quarks and leptons 



is increased from three to four or more [p34H237| . |l 58| . [238l , |239fl . In particular, Harari and 

Leurer |[238|| claim to have an optimal parametrisation for general generation number. Let 

us here merely quote ||158|| an example of a parametrisation for Vq KM in terms of six mixing 

angles and three CP-violating phases: 
/ 



(4) 



CKM 



Ci S1C3 

-sic 2 C1C2C3 + s 2 s 3 c 6 e z 



S1S3C5 
C1C2S3S5 - s 2 c 3 s 5 c 6 e 



iSi 



S1S3S5 

C1C2S3S5 - s 2 c 3 s 5 c 6 e l 
-s 2 s 5 s 6 e i{5l+S3) 



— S1S2C4 C1S2C3C4 — c 2 s 3 c 4 CQe l ^ x C1S2S3C4C5 + c 2 c 3 c 4 c^CQe lSl Cis 2 s 3 c 4 c^ + c 2 c 3 c 4 c^CQe 



-s 3 s 4 s 6 e 



~c 2 c i s 5 s 6 e i ^ +5 ^ 
+c 3 s A c 5 s 6 e lS2 
+s 4 s 5 c 6 e l(52+S[i) 



— S1S2S4 C1S2C3S4 — c^^s^CQe^ 1 C4S2S3S4C5 + c 2 c 3 s 4 c§e 



i8i 



+s 3 c 4 s 6 e 



iS-2 



-c 2 S4S 5 s 6 e 



i(Si+S 3 ) 



-c 3 c 4 c 5 s 6 e 



i&2 



-c A s 5 c 6 e 



HS2+S3) 



+c 2 c 4 s 5 s 6 e l{Sl+5 ^ 
+c 3 s 4 c 5 s 6 e lS2 
-s 4 s 5 c 6 e i(52+s ^ 

C1S2S3S4C5 + C 2 C 3 S 4 C^ 

+c 2 s 4 s 5 see^ 1+s ^ 

-c 3 c 4 c 5 s 6 e l52 
+c 4 s 5 c 6 e i(52+5a) 
(83) 



The large number of phases and mixing angles makes a thorough analysis of CP violation 
in the model with four chiral generations impractical. However, if one considers the case of 
an additional isosinglet quark, then the discussion of CP violation becomes tractable, and 
this case has all of the essential features of the more general case. An extremely extensive 
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and detailed analysis of CP violation in models with a Q = — 1/3 isosinglet quark has been 
given in a series of papers by Silverman ||166| , |170| , |161| , |162|| ; the Q = 2/3 case was discussed 



briefly by Barger et al. ||154j| . Using the notation of Barger et al. (see the last subsection of 



Section III), one can write the unitarity relations of the 4x4 matrix as 

KiVuj + V*V cj + V£V tj + V^oj = (84) 



or 



V* d V jd + V*V js + V* b V jb + V* D VjD = k 3 (85) 

For i ^ j, these can be expressed as closure of a unitarity quadrangle. The first three terms 
are the three sides of the unitarity triangle (which would be closed if the 3x3 submatrix were 
unitary). Note that in the first relation, the fourth side of the quadrangle is V^Voj = — 
which is the coefficient governing the size of the flavor- changing neutral current interactions. 
There are two major consequences for CP- violation in the B sector. The CP violating angles 
a and (3, which are directly measureable in B d decays, no longer have the same values as 
they do in the standard model; and tree-level Z-mediated graphs give a contribution to 
B d — overlineB d mixing. The contribution of z^b to the unitarity quadrangle has a fixed 
magnitude (since the quadrangle must close), but has any phase. Barger et al. give the CP 
asymmetry for B d — > ipK s as a function of two parameters: the phase of z bd and \z bd / {y td V^ b \. 
For some values of the parameters, the CP asymmetry is the same as the standard model; 
for some values it is different, and can even have the opposite sign. 



The most extensive analysis of the model is the recent paper of Silverman ||166|| . He 
first finds the presently allowed ranges for the unitarity triangle angles, the mixing and 
asymmetry in B s — B s mixing. Then, the model with an isosinglet Q = —1/3 is presented. 
He analyzes all of the constraints in this mdoel, including experiments to determine the 
CKM submatrix elements, |e|, Kl — > ~, B d — B d mixing, B s — B s mixing, B — > fi + fi~X, 
Rb in Z — > bb, and the recent event in K + — > 7i + uV. The results are presented as plots 
(one for the standard model and one for the isosinglet model) in the (p, rj) plane, in the 



(sin(2a), sin(2/3)) plane, and in the (x s ,sm , y) plane. In each case, he shows the current 
bounds (with la, 2a, 3a contours), and then shows projected results from the upcoming B- 
factories. There are several interesting features of the isosinglet model case: the presently 
observed CP violation in e can come entirely from new phases, and the value of sin 7 can 
take on any value whatsoever. The range of Ab s is just as large for its negative values as for 
its positive values. It is clear from the plots how the upcoming B-factories will drastically 
narrow the allowed parameter space, can could easily rule out standard model CP-violation. 



E. CP in the Aspon Model. 



Here we shall concentrate specifically on the Aspon Model [ |61| , p^ . p40| - |243| . |245| , p46| , p!8S 



and its requirement of additional quarks to solve strong CP. The first generation of the 
standard model contains quarks with the following (T 3 ,Y) values: 

1 U d L ;(o,l) = d L ; (86) 



2' 6J ' V ' 3 

i^) = ^ ; (°--|) = fi - < 87) 

We introduce |)T| a U{l) new symmetry and assign Q new = to all of the above quark states 



and to the leptons, although the latter do not play a significant role in solving strong CP. 
The second and third families have parallel assignments under the same U(l) new . 

In the model there is also a real representration of exotic heavy quarks corresponding 
to a complex representation C and its conjugate C. In C the exotic heavy quarks have 
quantum numbers exactly like some of the usual quarks; for example, in C there may be 
one doublet 

U L (88) 

These have charges Q new = +h. In representation C we shall then have 

1 ^ ~,C.f 1 K\_ TT C 



2' 6 )- D t\-2>-i>)= V i < 89 > 
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These have Q new = —h. 

The Higgs sector has one complex doublet 

1 1 







2' 2 



i Qnew 0, 



(90) 



and two complex singlets 



Xi,2(0,0),Q neu) = +h. 



(91) 



The gauge group is SU(3)c x SU(2) L x C/(l)y x £/(l) ne u>, where we gauge £/(l) new to avoid 
an unwanted Goldstone boson when it is spontaneously broken. 

In breaking the symmetry, we give a real vacuum expectation value to <fi and complex 
VEVs to xi,2 with a nonvanishing relative phase. 

The lagrangian contains bare mass terms M(U^Ul + D^Dl) for the extra quarks. The 
allowed Yukawa coupling include u l L u J L <f), d l L d J L 4> for the families and U^u l L Xa, D^d^Xa cou- 
pling light quarks to C heavy quarks (a = 1,2 and i = 1,2,3). Because the families have 
no couplings to C exotics, the quark mass matrix determinant arising from spontaneous 
symmetry breaking is real at lowest order; it has the required texture. 

We do not allow terms in the Higgs potential which explicitly break U(l) new . Disallowed 
terms include </>0x 2 > X 2 > X 3 an d X 4 - If an y °f these terms are present, U(l) new is explicitly 
broken and the model can have = at tree level only in very special cases where e.g. 
we choose particular representations of a grand unified group such that the quark matrix is 
real. 

Without explicit breaking of the U(l) neW) there is correct texture at tree level; the mass 
matrix has the tree-level texture (F = family) 

/ 



{FCC) 



V 



c 



(92) 



real complex 
complex real 
real 

Thus Qqcd — at tree level. If we assume CP symmetry of the Lagrangian then 
©qfd = also. In this case = at tree level and will be nonzero by a small amount 
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through radiative corrections; this can be consistent with experiment if the Yukawa couplings 
are within a certain window, as will be shown below. 

Because it is anomaly-free, we may gauge the U(l) new symmetry and the Higgs mecha- 
nism will lead to a massive gauge boson, the aspon, which couples only to the exotic quarks 
and indirectly via nondiagonal mixing to quarks and leptons. 

The Yukawa interactions of the model are given by 



L Y = qLm d d R 



V2 



+ q L m u u R 



$ 

v 



V2 



+ ti*q L Q R x a + h.c. (93) 



where v/y/2 is defined as the VEV of 0° and $ as (0°, — 4>~) T ■ The generation indices are 
implicit. Usual quarks and leptons acquire their masses through spontaneous symmetry 
breaking (SSB) induced by the VEV of the doublet Higgs scalar. The new quarks acquire 
their mass through a gauge invariant mass term of the form MQ L Q R . Hence, U and D 
quarks are degenerate in mass, m^, m u , m e , v, h 1 ' 2 , and M are real by the assumption of CP 
invariance. The VEVs of xi an d X2 are chosen to be 



< Xi >= ^ Kie 5 <X2>= -J=K2 



(94) 



Hence CP is broken spontaneously. [CP can be broken softly by i(x*X2 — X2X1)] 

The up- and down- quark mass matrices linking the right-handed sector to the left-handed 
sector are of the form 



M, 



up 



m„ F 

M 



M, 



down 



m d F 

M 



(95) 



where 



h 1 < Xi > +h 2 <X2> 



(96) 



The Kobayashi-Maskawa (KM) matrices will be generalized to 4 x 4. From the constraint 
|Kd| 2 + |K s | 2 + |K6| 2 = 0.9979±0.0021, we find \F X \/M and \F 2 \/M to be less than lO" 2 and 



10 , respectively. Although F is a complex column matrix, the determinants of M up and 
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Mdovm are real. All entries become complex and, therefore a nonvanishing value of arises 
through radiative corrections. The calculation of at one loop level will be done below. 

First, we discuss how the flavor-changing neutral currents (FCNC) in the presence of 
new quarks are suppressed at the tree level. After introducing the new vector-like quark 
doublet, we find that there are FCNC's induced by Z coupling because of the mismatch 
of the new and usual quarks in the right-handed sector. Therefore, the flavor-changing Z 
couplings are induced by the terms 

1\ 92 



FCNC 



-Dr^^DrZ 11 + (Ur contributions), 



(97) 



2 J cosOw 

where the factor — \ is the isospin of D R and g 2 is the SU(2) gauge coupling constant. Let 
us consider the down sector first. Without losing any generaity, we assume the down-quark 
mass matrix is in the partially diagonalized form 



M, 



down 



m d Fx 

m s F 2 

m b F 3 

M 



(98) 



This mass matrix can be diagonalized by a biunitary transformation, K^M^ ou , n K^. The 
transformation matrices are given in Ref. p2]| . Thus Eq.fl9?D can be written in terms of 
mass eigenstates d' 1 as 



L F z CNC (down) = P^kludlZ, 



(99) 



for i ^ j, where 



Pi 
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-(K R y 4i (K R ) 4j 



2 J cosOw 
1\ 92 rn di m dj 

2 J cos6 w M 2 1 j 



(100) 



where Xi = Fi/M. Therefore, the FCNC induced by Z coupling is highly suppressed by the 
small mass ratio of usual to new quarks. It is because the mixings of right-handed quarks 
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require a helicity flip of the usual quarks. For example, ~ 544.8 x 10 _8 xiX2 < 10~ 10 for 
M = lOOGeV, while the experimental limits on FCNC's require only that /3y < 10~ 6 . 
FCNC can also be induced by aspon (A) couplings which are given by 



L F A CNC {down) = aijD'l^dlA" 



(101) 



for i ^ j, where 



a; 



9A%i% J 



(102) 



Therefore, FCNC's induced by A will be important if A is not too heavy compared to Z. 
Consider the K° — K° mixing matrix element M 12 . e(M 12 ) is expected to be dominated by 
standard 21^-exchange box diagrams, while Im(Mi2) receives its largest contribution from 
the A exchange shown in Figure 16. 




FIG. 16. Contributions to Im(AMi2) by aspon exchange (Circles mean mixings.) 



We obtain 



Im{M_ 



12 J 



ffrriK 1 

6 K 2 



Im(xiX* 2 ) 2 



(103) 



where k 2 = k\ + k, 2 ,. The color factor has been taken into account in Eq.( 103 ). Im(M 12 ) 
receives contributions from the new 2W / -exchange box diagrams shown in Figure 17, 
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(a) 



■i o°- 



u 



w 



w 



-o- 



u 

(b) 

FIG. 17. Contributions to Im(AMi2) by new quarks and two-W box diagrams. (Circles mean 
mixings.) 

but these contributions are negligible. We will consider the CP-violating parameters, such 
as Im(Mi2) / AMk in more detail below. 

Next we consider the up sector for completeness. We can also choose states such that 
M up in Eq.(|95|) is replaced by 

^ m u 



M, 



up 



m c F 2 
m t F 3 
M 



(104) 



with 



F = C F 



(105) 



where C is the real standard 3x3 KM matrix. The transformation matrices J l and Jr that 
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diagonalize M up in Eq. (|104| ) can be related to and K/j by changing Xi into Xi(= Fi/M) 
and rrid, m s and mb into m„, m c and m t . The generalized 4x4 KM matrix is given by 



V 



(4) 
KM 



Jt 



C 
1 



K, 



(106) 



Before proceeding, we note that flavor-changing Z-coupling can be induced by the one-loop 
diagram shown in Figure 18. 




FIG. 18. Flavor-changing Z coupling induced by new quarks at one loop. 



By naive dimensional arguments, the effective coupling constant is 



1\ g 2 hjhj M_ 
2/ cos9 w lQn 2 \m Xj 



(107) 



Using /ii,2,3 — 0.01 and (M/m x ) 2 ~ 0.1, we conservatively estimate /3y to be less than 10 -7 . 
Therefore, we expect these FCNC's to be smaller than those in the standard model. 

Consider now the one- loop corrections to 0. Although the mass matrices in Eq.([55|) are 
complex, their determinants are real. Therefore, defined in Eq.(|8"0D is zero at tree level. 
will be nonzero when the mass matrices receive radiative corrections. For example, the 
contributions to 9 from the up sector are given by 

e(up) = Arg[det(M up + 5M up )} 



ImTrln[M up (l + M^SM^)} 
ImTr{M-^M up ), 



(108) 
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where we have used the fact that M up is real and that the corresponding radiative corrections 
5M up are small. The last line in Eq.( |108| ) is valid at one- loop order. Defining the one- loop 
corrections 5M up by 



6M, 



1 5m„ 5F ^ 



up 



5m Uu 8M 



(109) 



and combining with Eq.(|108|), we obtain 



G(up) = ImTr(ra u l 5m u - m^FM^Smuu + M~ X 8M). 



(110) 



Notice that SF will not contribute to at one-loop order because of the structure of M up . 
The expression for Q(down) is strictly analogous to Eq. flllOj ). 

By studying all possible one-loop diagrams |62lJ246|1 we find that the only contribution to 
G comes through 5m u in Eq. (|110|) and, in particular, from the diagram depicted in Figure 
19. The imaginary part gives a contribution 

a=21=3 



G(up) = — \= —, — ~ttt Y h?Im[x*]\ a ^- 



(111) 



(<t>> 



\ / 



/cj> 



u 



R 



R 



Ui 



FIG. 19. One loop diagram of which the imaginary part contributes to O. 
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which gives the estimate (for convenience we shall write x to denote x = \xi\, the modu- 
lus, and taken to be generation independent since the limits on x are not sensitive to the 
generation considered) 

Xx 2 

Here X a is the coefficient of the quartic interaction between the two types of Higgs |0| 2 |x a | 2 > 
and A with no subscript is an average value. (Actually there are three independent A 
corresponding to indices 11, 12+21, 22 but our estimates will not distinguish these). 

As mentioned earlier, the neutron electric dipole moment d n has been calculated [ [213|J214| 
in terms of G long ago with the result that 

d n ~ 10~ 15 6e.cm. (113) 

and so we know from d n < 10~ 25 e.cm. empirically that B < 10 -10 , from which it follows 
by Eq.([TT2|) that Aa; 2 is less than 1(T 8 . 



In the kaon system, the CP violation parameter |e^| is given p2"| , |245| by 

Using (Am^/m^) = 7.0 x 10~ 15 , fx = 0.16GeV gives the relationship between x 2 and the 
U(l)new breaking scale 

k/x 2 = 2.9 x 10 7 GeV (115) 

Thus, if we insist that the U(l) new is broken above the electroweak breaking scale (~ 
250GeV) then 

x 2 > 10^ 5 (116) 
From Eq. ( |112| ), this means that A < 10~ 3 . 



In [ f243[ | , it was argued plausibly that A > 10 5 on the basis of naturalness; this would 
imply that > 10 -12 and hence d n > 10 _27 e. cm. 
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But we can find a yet more solid and conservative lower bound on d n . It follows from 
the fact that the |</>| 2 |x«| 2 interaction receives a one- loop correction from the quark loop 
box diagram where three sides are the top quark and the fourth is the heavy [/-quark (see 
Figure 20). 



X 



U R 



"R 



X 




4> ♦ 

FIG. 20. One loop |0| 2 |x| 2 counterterm for the coupling A. 



The full A is given by 



A = ^tree{bare) + Ai_z oop (including counterterm) + higher loops 



117) 



and the one- loop finite contribution, for the dominant diagram, Figure (20), neglecting 
quark masses and taking /if, g t as the respective Yukawa copuplings to Xa and <f> of the third 
generation by 



^i-ioop ^ J 7^Xil^3 1 2 \9t \ 2 jj + counterterm 
which imples that the lowest value for A (without accidental cancellations) is: 

.,2 



(118) 



A > 



x- 



16tt 2 



(119) 



Combining Eqs.( 112 ) and (119) then gives the estimate for of 



e > 



X 

16vr 4 
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(120) 



which implies that x 2 < 10 3 (incidentally in full agreement with ||243|| ) and that 10 10 > 

> io- 14 . 



From Eq.( |113| ), this then gives a lower limit on the neutron electric dipole moment of 

d n > 10~ 29 e.cm. (121) 

This is more than two orders of magnitude greater than the prediction of the KM mechanism 
and thus provides another distinguishing feature of spontaneous CP violation. 

Before addressing other consequences of the Aspon Model, let us point out that the 
production of the predicted heavy quarks and the aspon in a hadron collider is discussed in 
Ref. pq| . 



A promising approach to finding new quark flavors is through searching for heavy quark 
bound states. This fails for the top quark since the single t quark decay to bW is more 
rapid than the formation time of toponia. However, the Q of the Aspon Model have a very 
much longer lifetime because of the small mixing with the three light generations. Such 
production of the exotic quarkonia is dicussed in . 



The experimental information regarding CP violation still comes primarily from the 
neutral kaon system and is inadequate to determine whether the KM mechanism is the 
correct underpinning of CP violation. In dedicated B studies, with more than 10 8 samples 
of B°(B°) decay, it will be possible ||247| - |25l| to test this assumption stringently by measuring 



the angles of the well-known unitarity triangle whose sides correspond to the complex terms 
of Eq. (|78|) . If CP is spontaneously broken, as in the Aspon Model, the outcome of these 
measurements will be different from the predictions of the standard model. 

The three angles of the unitarity trangle are conventionally defined as a, (3, and 7 be- 
tween the first and second, second and third, third and first sides in Eq.(|78|) respectively. 
These angles can be separately measured for the standard model by the time-dependent CP 
asymmetry 

_ T(B°(t) -> f)-T(B° -> f) 

af[t) ~ r(B°(t)^f) + r(B°-^f) (122) 
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where the final state / is a CP eigenstate. We define q,p in B° — B° mixing by the mass 
eigenstates B 12 : 

\B xfl >= p\B° > ±q\B° > (123) 

and similarly for K\ 2 in the kaon system. Also, A, A are the decay amplitudes 

A, A =< f\H\B°,B° > (124) 

. Let us consider the specific cases of / = ti~ + it~ , ipKs from Bd decay and / = pKs from B s 
decay. We define A(/) by 

A(7r + = (-1 (4) > ( 125 ) 

\p)b a \ A ) Bd ^- 



K*Ks) = [±\ (4l (A , (126) 



and 



KPK S )=(1) 4 (?) . (127) 



\P/S fl \ A /B.-pK \P/ K 

The complex conjugate appears in Eq. (|127|) because 5° — > K° while 5° — ► if . If to a 



sufficiently good approximation \q/p\ = 1 and [v4./v4.| = 1, as we shall show for the Asp on 
Model below, then A(/) is related to the CP asymmetry through the B\ —B 2 mass difference 
AM by 

a f (t) = -lm\(f)sm(5Mt). (128) 

In the standard model the angles of the unitarity triangle are related to the A(/) by: 

sin2a = ImA(7r + 7r~); sin2(3 = lmX(ipK s ); sin27 = lm\(pK s ). (129) 

Such relations are no longer valid in the Aspon Model because Im(q/p)B d has a major 
contribution from aspon exchange and Im(q/p)K is dominated by aspon exchange. 
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To evaluate the CP asymmetries in B decays for the Aspon Model one needs to evaluate 
the different factors in the A(/) given in Eqs.( |125| ) - ( |127| ) above. More precisely we need, 
from Eq. ( |128| ), the imaginary part of the A(/). The Aspon Model adds new Feynman dia- 
grams involving aspon exchange to those already present in the standard model. Because 
CP is only spontaneously broken, the W-exchange amplitudes are predominantly real amd 
have very small phases while the aspon exchange has a much smaller magnitude but an un- 
predicted arbitrary phase. As a result, the |ImA(/)| appearing in Eq.( |128| ) are of order 0.002 
or less, compared to the standard model expectation that |ImA(/)| be of order of, although 
less than, unity. Thus CP asymmetries in B decays are predicted to be correspondingly 
smaller than in the standard model. The technical details can be found in Ref. p42j| . 

It is found, setting M A = 300 GeV, that 

|ImA(7r + 7r-)| < 1 x 10~ 5 , (130) 
\Im.X(ipK s )\ < 2 x 10~ 3 , (131) 

\lm\{pK s )\ < 2 x 10" 3 . (132) 

The resulting asymmetries a/(t) are much smaller than those predicted by the standard 
model where these imaginary parts are all of order unity. 

It is also interesting to observe from Eqs. (|125|) - ( |127f ) that 



\{jjK s )\{ P K s ) = fq\ (A\ 

A(7T + 7T-) " [p) Bs \A) Bd ^ K 

= - 

\Pj B s \ A J B S ->D+D7 

= X(DfD;). (133) 

In the Aspon Model, where the A(/) have unit moduli and |ImA(/)| 1, this relation 
implies a linear relation for the imaginary parts: 

lm\(ijK s ) + lmX(pK s ) - ImA(7r + 7r-) - lmX(D+Dj) = 0, (134) 
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which provides an additional test of the Aspon Model. 

In conclusion, our result is that, if the Aspon Model is correct, CP asymmetries in B 
decays would be much smaller than predicted by the standard model and the relation (|134|) 
would be satisfied. 

One may ask | |243| ] whether the present experimental situation of B decay is compatible 



with the Aspon Model? 

In the model, the quark mixing matrix is a complex unitary 4x4 matrix C av . Mixings 
of the conventional six quarks with one another are specified by the 3x3 matrix CV,-, whose 
indices run from 1 to 3. CV,- is neither real nor unitary because of \~ induced mixings to 
the undiscovered quark doublet. However, these terms are ~ x 2 . Thus, since x 2 < 1CT 3 , 
Cij is, to a precision of at least 0.1%, a real orthogonal matrix. It is a generalized Cabibbo 
matrix rather than a Kobayashi-Makawa matrix. This is unfortunate for the search for CP 
violation in the beauty sector, but has other observable consequences. 

In the standard model, the Kobayashi-Maskawa matrix V is complex and unitary. The 
sides of the unitarity trangle are unity and 



Ri, 



VudVub 


, Rt — 


V td Vtb 


VcdVcb 


VcdV cb 





(135) 



The value of has been measured. According to ||251]| , Rj, = 0.35 ± 0.09. The value of R t 
cannot be extracted from experimental data alone. Appeal must be made to a theoretical 
evaluation of the neutral B-meson mass difference using the standard model. The analysis in 
| 251|| yields Rt = 0.99 ± 0.22. These results suggest a rather large value of the CP-violating 



angle /3, namely, 0.34 < (3 < 0.75. 

In the Aspon Model, the matrix is orthogonal up to terms of order x 2 arising from 
mixings with unobserved quarks. Thus, we anticipate no readily observable manifestations 
of CP violation in the beauty sector. Furthermore, the unitarity triangle must degenerate 
into a straight line: \Rb±Rt\ = 1. In this case, we cannot appeal to a theoretical calculation 
of the neutral B-meson mass difference since it depends on unknown parameters. On the 
other hand, the matrix with neglect of terms ~ x 2 involves only three parameters. ^From 
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data at hand, in this context, we obtain R t = 1 — p in the Wolfenstein [|252|| parametrization, 
whence R t = 0.637 ± 0.09. We note that present data yields R b + R t = 0.99 ± 0.13. This 
result is compatible with an approximately orthogonal mixing matrix and hence with the 
Aspon Model. 

As a final CP violation parameter in the Aspon Model, we shall discuss the value pre- 
dicted for Re^^ - a measure of direct CP violation in K decay. To evaluate Re^^j 
requires the study of several Feynman diagrams j253H258H, and their comparison to the 
standard model. Recall that the most recent evaluations completed at CERN (NA31) 
1 259|| and FNAL (E371) |260|j26T| give results Re (4) = ( 23 ± 3 ' 6 ± 5 ' 4 ) x 10 ~ 4 and 
Ref—^ = (7.4 ± 5.2 ± 2.9) x 10 -4 respectively, where the first error is statistical and the 
second is systematic. These results are consistent within two standard deviations; the error 
is expected to be reduced to 1 x 10~ 4 in foreseeable future experiment sf]. 

A detailed analysis of e'/e in the Aspon model was carried out by Frampton and Harada 



245 , |246|| . They showed that penguin diagrams involving the additional quarks give the 



dominant contribution. The outcome of these considerations is that Rel ~ J is not larger 
than 10~ 5 in the Aspon Model. While it does not vanish exactly, it does correspond closely 
to a superweak model prediction |264 . 

This discussion of the Aspon Model is presented as a motivation for additional quarks 
beyond the six discovered flavors. Because of the smallness of x 2 which characterizes the 
mixing of the additional quarks with the known ones, the new quarks have a long lifetime. 
This is important in their experimental detection, as discussed the next Section. 



6 A new result from KTeV |244) | at Fermilab gives Re (^j = (28 ±4.1) x 1CT 4 ; it will be interesting 
to see whether CERN confirms this result. 
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VII. EXPERIMENTAL SEARCHES. 



A. Search for long-lived quarks 



1. Present Searches 



The search for long-lived quarks is an ongoing process at the Fermilab Tevatron. The 



first search at the Tevatron was made by the DO collaboration [265] which looked for signals 



b' — > by, where b' is a charge -1/3 quark, setting a limit my > Mz + m^. The second and 



most recent search was made by the CDF collaboration ||266|| which looked for a displaced 
vertex for Z — > e + e~ coming from b' — > bZ resulting in my > 148 GeV (for cr = 1 cm.). The 
latter search will be described below. Eventually, such a search will be carried out at the 
LHC which has a much greater CM. energy, with a much larger production cross section. In 
this section, we will concentrate on the current limits coming out of two operating facilities: 
LEP2 and Fermilab. In the next section, we will discuss the prospects for future searches, 
both at the Fermilab Tevatron and at the LHC. We will only briefly discuss prospects for 
such a search at facilities which are under discussion, such as the NLC, etc. Also in this 
section, we will focus primarily on hadron colliders such as the Tevatron since these are the 
machines which can explore the mass range that was discussed earlier. 

In the search for a new particle, there are two principal activities to attend to: how 
to produce the particle and how to detect it. For a particle which is somewhat "exotic", 
such as supersymmetric particles, the production process would be highly model-dependent. 
Fortunately, for a heavy quark, this is rather standard: it proceeds through the qq and gg 
channels. For the range of heavy quark masses considered in this Report, the qq process via 
the electroweak channels W, 7, Z is completely negligible compared with the QCD process 
with gluons. The production cross section, at the Tevatron and at the LHC, for the top 
quark as a function of its mass has been computed up to the next-to-leading order in QCD 
|267|| - |[272| . (The qq channel is dominant at the Tevatron while the gg channel is dominant 



at the LHC.) This can be directly applied to the present search for long-lived quarks whose 
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production mechanism should be similar to that of the top quark. The production cross 
sections as a function of the heavy quark mass are shown in Figs. 21 and 22 for the Tevatron 
at y/s = 1.8 TeV and 2 TeV respectively, and in Fig. 23 for the LHC at y/s = 10, 14 TeV. 
Here "m t " will stand for a generic heavy quark mass, for both the Tevatron and the LHC. 
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FIG. 21. Physical cross section for — > ttX at ^/s = 1.8 TeV as a function of the top mass. 
This cross section applies to a heavy quark Q as well with t changed to Q. The two data points 
are from CDF and DO respectively for the top quark. 
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FIG. 22. Physical cross section for pp — > ttX at yfs = 2.0 TeV as a function of the top mass. 
This cross section applies to a heavy quark Q as well with t replaced by Q. 
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FIG. 23. Cross section for pp — > ttX at yfs = 10 and 14 TeV as a function of the top mass. 
Notice that p has been mislabeled as p in the figure. The same prediction applies to a heavy quark 
Q. 

As can be seen above, the predicted cross section at the LHC, for a given heavy quark 
mass, exceeds that at the Tevatron by more than two orders of magnitude, which will 
facilitate the search for such an object. 

The next task is to define the detection capability of various detectors. Since the lat- 
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est constraint on long-lived quarks come from CDF [ 266|| , we shall use it as a prototype 



of detectors dedicated to such a purpose. Other detectors such as DO are very similar in 
layout. Needless to say, the CDF detector is a complicated, multipurpose one whose spec- 
ifications can be found in |[273|| . A generic detector for hadron colliders generally consists 
of a (silicon) vertex detector immediately surrounding the beam pipe for high precision de- 
termination of the location of the tracks. Next comes a central tracking chamber which 
measures charged tracks and momenta of charged particles. Surrounding these two units 
are generally hadron calorimeters which measure the energy deposited by hadrons. Next 
comes the muon chambers which detect the location of the particles which penetrate the 
calorimeters. As described below, the first two parts (vertex detector and central tracking 
chamber) were used to search for displaced vertices coming from the decay of a long-lived 
particle. For stable or very long lived particles, the muon chambers are used in conjunction 
with the ionization energy loss in the tracking chambers to make such a search. 

The current search at CDF can be divided into two categories: the search for those 
quarks whose decay lengths, / = 7/3cr with r being the proper decay time, is 1) less than 1 
meter, and 2) greater than 1 meter. 

Let us first concentrate on the first category (/ < 1 meter) | |266[ | . The parts of the detector 
which are relevant here consist of two components: a silicon vertex detector immediately 
surrounding the beam pipe for precision tracking and a central tracking chamber embedded 
in a 1.4 T solenoid magnetic field which measures the momenta and trajectories of charged 
particles. In the search for a new particle, a crucial task would be the identification of a 
characteristic signature which would distinguish it from background. In the present case, 
that characteristic signature is the decay Z — > e + e~ with the e + e~ vertex displaced from 
the pp interaction point. This Z boson could come from the decay of a charge -1/3 quark 
(denoted by b' in Ref. [ |266|| and by D in Ref. ||188|| in the process D — > b+Z with a subsequent 



decay Z — > e + e~ . In this case, D would be the long-lived parent of the Z boson. The search 
at CDF concentrated on events containing an electron-positron pair whose invariant mass is 
consistent with the Z mass and whose vertex is displaced from the pp interaction point. The 
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data used was from the 1993-1995 Tevatron run with an integrated luminosity of 90 pb^ 1 of 
pp collisions at y/s = 1.8 TeV. 

In the search for a long-lived parent of the Z, the CDF collaboration focused on the 
measurement of L xy which is the distance in the transverse (r — 0) plane between the pp in- 
teraction point and the e + e~ vertex. Notice that L xy = ^f3 xy CT with (3 xy being the transverse 
component of the parent particle didived by c. As defined by the CDF collaboration, L xy can 
be either positive or negative. For prompt Z's coming from the SM process qq — > Z — > e + e~, 
one would expect L xy m because of the short lifetime of the Z. The L xy distribution with 
appropriate cuts taken into account is shown in Fig. 24 below 




FIG. 24. The L xy distribution of the Z's after applying all cuts. The data are represented by 
the circles. The histogram is the expected L xy distribution for prompt Z's based on the measured 
L xy uncertainty in the event sample. The inset shows the distribution after the 2 jet requirement 
is applied. The vertical dashed lines separate the prompt and non-prompt regions. 
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As emphasized by [ ^66] , this distribution is consistent with that for prompt Z's where one 



would expect less than one event for \L xy \ < 0.1 cm. ||266|| also pointed out that the number 
of events with L xy significantly less than zero is an effective measure of the background. The 
CDF collaboration observed one event for L xy > 0.1 cm and 3 events for L xy < —0.1 cm. 
As stated, there is no evidence for a long-lived parent of the Z. This is shown in Fig. 25 
where the constraint is expressed in terms of the 95 level upper limit on the product of the 
production cross section for the long-lived parent, ax, its branching ratio, Br(X — > Z), the 
branching ratio, Br(Z — > e + e~), and the e + e~ acceptance for pseudorapidity \r)\ < 1. 
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FIG. 25. The 95% confidence level upper cross section limit for a. Br times the acceptance for 
an electron-positron pair to be within the detector as a function of fixed \ xy = "if3 xy CT. Cross 
sections above the curve have been excluded at the 95% confidence level. The inset shows the 
exclusion curve and the theoretical prediction for a b' quark of mass 110 GeV as a function of its 
lifetime, assuming 100% decay into bZ. 
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The above discussions and figures deal with limits on the production of a single parent 
with its subsequent decay into a Z. We are, however, most interested in the detection of a 
long-lived quark. As we have mentioned earlier, this long-lived quark would be produced in 
pair. The CDF search for a long-lived D (or b') which is pair-produced can be summarized 
as follows. The kind of events that are searched for would include, besides the e + e~ pair 
coming from the Z, two or more jets. For instance, this could come from the reaction: 
qq — > DD — > bZbZ — > be + e~bqq. The L xy distribution is shown in the inset of Fig. 25. 
There one would expect less than one event for L xy < 0.01 cm. The CDF collaboration 
found one event. This was then translated into a cross section limit as a function of cr, 
where r is the lifetime of D. Assuming that Br(D — > bZ) = 100%, the exclusion curve 
for the DD production cross section as a function of cr is shown as an inset of Fig. 25 
for a D quark of mass 110 GeV. The theoretical prediction for the cross section for such a 
mass is shown as a horizontal line. Clearly, this is ruled out for a wide range of lifetimes. 
For other masses, the exclusion curves for the cross section are not shown but are instead 
translated into exclusion regions in the mass-lifetime plane. (This is because the production 
cross section can be calculated in QCD as a function of the D mass as mentioned above.) 
The plot shown in Fig. 26 assumes the above branching ratio. 
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FIG. 26. The hatched areas in this plot represent the 95% confidence-level regions of b' mass 
and lifetime that have been excluded. For cr = 1 cm, CDF excluded up to a mass of 148 GeV. 

In Fig. 26, three forbidden regions are presented: the LEP, the -DO, and the CDF 
constraints. The most stringent constraint comes, of course, from the CDF results. As can 
be seen from Fig. 26, for every cr, there is a range of forbidden masses represented by the 
shaded region. The largest forbidden range is for cr = 1 cm corresponding to a lifetime 
r m 3.3 x 10~ n sec . This rules out the mass of the D quark up to 148 GeV. For smaller or 
larger cr, we can see that the lower bounds on the D mass become somewhat smaller than 
148 GeV . If the D quark, happens to have a mass larger that 148 GeV, it could escape 
detection for a large range of lifetimes as can be seen in Fig. 26. The mere fact that there 
exists unexplored regions of the detector, as shown in the unshaded areas of Fig. 26, is 
reason to believe that there are plenty of opportunities for future searches. 

What does this result tell us about the long-lived quarks with the mass range that we 
have discussed in Section IV? First, as we have mentioned earlier, we are especially concerned 
with the mass of the heavy quarks larger than 150 GeV. From Fig. 26, one can see that the 
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CDF search based on the decay mode D — > hZ does not set any constraint on long-lived 
quarks with a mass greater than 150 GeV. In other words, if this long-lived quark exists 
and if it decays at a distance L xy > 0.01 cm, it has yet to be discovered. This statement 
is, of course, based on the assumption that the branching ratio for D — > bZ is 100%. As 
discussed in Section IV and in Ref. [|188| , for m D < m t , there is another possible decay mode 



for D, namely D — > (c, u)W (the c quark channel in any reasonable scenario dominates over 
the u quark channel). Whether or not D — > bZ dominates over D — > cW will depend on 
a particular model for the mixing element |Vd c |. As discussed in Ref. [ p. 8 8(1 , even with a 
very naive assumption |Vb c | ~ x 3//2 , with x being the mixing parameter between the third 
generation and the heavy quark, D — > bZ can dominate over D — > cW for a certain range 
of mass and mixing x. For |Vd c | ~ x<2 ( or less), the bZ mode will almost always be the 
dominant one. Since this is a rather model-dependent statement, one should, in principle, 
look for both modes when vhd < m t . Unfortunately, the mode D — > cW would be rather 
difficult to detect. We shall come back to this issue and others in the discussion of future 
searches. 

The next question concerns the limits on a charged "stable" or very long-lived massive 
particle. Such a search is being carried out a CDF. Basically, this search focuses on decay 
lengths larger than 1 m, i.e. larger than the radius of the Central Tracking Chamber. 
As of this writing, preliminary results have only appeared in conference talks, |113| , |114| . 



Therefore, what is described below will be considered preliminary. A stable massive quark 
moving at a low velocity will leave an ionization track in the tracking chamber because 
the energy loss dE/dx oc l/(3 2 (the Bethe-Bloch equation) and a low (3 would imply a 
large energy loss. Furthermore, for such a stable massive quark to be detected, one would 
look for signals in the muon detector after mesons formed from this particular quark have 
traversed the calorimeters and reached the muon detector. To distinguish it from a muon, one 
would have to correlate this signal with the large energy deposited in the tracking chamber. 
The measurement of dE/dx as a function of (3^ = p/M, combined with the momentum 
measurement would allow for a determination of the mass of the particle. The mass limits 
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from CDF for such a very long-lived quark are approximately 200 GeV. 

In summary, the most stringent limits, so far. on long-lived quarks come from the CDF 



collaboration ||266|| . It excludes a long-lived, charge -1/3 quark of mass up to 148 GeV for 
a lifetime r w 3.3 x 10~ n sec (cr = 1 cm). For other values of lifetimes, the excluded 
mass ranges are weaker as can be seen in Fig. 26. This constraint was based on the search 
for a displaced vertex for the decay Z — > e + e~ which could come from the decay D — > bZ. 
Furthermore, all of these constraints come from the search for the decay of the D quark inside 
the Central Tracking Chamber. If the D quark lives long enough to enter the calorimeters 
and subsequently trigger a signal in the muon chamber, the constraint (which is preliminary) 
is much stronger: a D quark mass below approximately 200 GeV is excluded. In short, under 
what circumstances will a D quark escape detection so far? First, if its mass is above 148 
GeV (as referred to in Section IV) and if it decays inside the Central Tracking Chamber. If 
the mass is above ~ 200 GeV, D is no longer required to decay in the tracking chamber: it 
simply escapes detection regardless of where it decays. Most of the discussion in Section IV 
concerned these possibilities. 

What kind of improvements should be made in order to be able to search for these quarks 
heavier than 148 GeV which could either decay inside the Central Tracking Chamber or, 
if heavier than 200 GeV, could also travel through the calorimeters? What about the U 
quark? How could one detect it? What if the dominant decay mode of the D is D — > cWl 
These are the kinds of questions that one would like to address, at least qualitatively, in the 
next section. 



2. Future Searches 

The first kind of future searches would be based on present facilities such as the Tevatron. 
In particular, one might ask what kind of improvement one can make by exploiting the 
present CDF RunI data with up to 120 pb~ l . One can then ask what light Runll with a 
large improvement in luminosity and detector might shed on the search for long-lived quarks. 
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The above search has focused on the discovery of displaced vertices for decay lengths 
greater than 100 /im, with the resulting constraints as described above. What if the decay 
length is less than 100 /xm? The appropriate kind of experiment would be a counting 
experiment which is not based on the search for displaced vertices |[274|| . How feasible this 



kind of experiment might be will probably depend on the improvements planned for Runll. 

These improvements include: a) a detector upgrade with, among several things, added 
layers of silicon; b) an increase by a factor of 20 in the luminosity. The radius of the silicon 
vertex detector is roughly of the order of 22.3 cm. Added silicon layers would increase that 
radius to about 28 cm [[274|| and consequently the tracking ability of the detector. 

For decay lengths between 100 fim and 1 m, one of the most important tasks would be 
to improve the tracking efficiency of the Central Tracking Chamber by adding,for instance, 
more silicon layers to the vertex detector. Decay lengths of a few tens of cm might be hard, 
although possible, to detect because of poor tracking efficiency in such a region. This would 
require new reconstruction algorithm. 

The search for very long-lived or "stable" quarks will also be improved by the detector 
upgrade and the increase in luminosity. 

One might ask what else could be done at CDF and DO in the next run beside those issues 
discussed above. In particular, one would like to know how feasible might the detection of a 
signal such as D — > cW be if it happens to be the dominant decay mode of the D. Needless to 
say, such a task would be much more daunting than the detection of D — > bZ. Nevertheless, 
a feasibility study would probably be extremely useful. As we have discussed at length in 
Section IV and in Ref. ||188|| , there is also the partner of the D, namely the charge 2/3 quark 



denoted by U, which should not be forgotten. If U is heavier than D- but not by much 
because of the p-parameter constraint- it will decay into D via U — > D + (/ + z/, g 2 / 3 g— 1/3), 
where particles inside the brackets denote light quarks or leptons or it can decay via U — > bW, 
depending on how degenerate U and D are and how large \Vub\ is. For the first mode, it 
was shown in Ref. [|188|] that U practically decays near the interaction point, with a decay 



length typically of the order of 10 5 /zm. The D will subsequently decay between 100/xm 
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and 1 m. It will be a challenge to be able to identify such a signal. For the second mode 
U — ► bW, one has to be able to distinguish it from a signal coming from top decay. It would 
be extremely hard, if not impossible, to be able to resolve the decay vertex to distinguish 
U from t. However, by comparing the predicted number of £'s with the observed ones, one 
might rule out the mode U — > bW with a U mass close to the top mass. 

Turning our attention to the upcoming experiments at the LHC, we would like to briefly 
describe the two main detectors which will be crucial to the search for long-lived quarks 
(if such a search would be carried out). They are the Compact Muon Solenoid (CMS) 
and a Toroidal LHC Apparatus (ATLAS) detectors [275]. The layout for both detectors is 
generically very similar to CDF and DO. 

In the search for long-lived quarks, the components which are crucial would be the central 
tracking system of CMS and the inner detector of ATLAS. The central tracking system of 
CMS consists of silicon pixels, silicon and gas microstrip detectors with high resolutions. 
(The resolution of the silicon pixels is about 11-17 fim while the outermost part of the 
detector, namely the gas microstrip detector, has a resolution of approximately 2mm!) The 
silicon pixels and silicon microstrips cover a radial region up to about 40 cm, a marked 
improvement over the CDF vertex detector. The microstrip gas chamber covers a radial 
region to approximately 1.18 m which is roughly similar to that covered by the Central 
Tracking Chamber of CDF. 

The inner detector of ATLAS consists of a Semi-Conductor Tracker (SCT): pixel de- 
tectors, silicon microstrips and GaAs detector, and Microstrip Gas Counters (MSGC). The 
pixel detectors have a spatial resolution of about 14 /im while the MSGC have a resolution 
of about 1.8 mm, very comparable to CMS. The SCT part covers a radial region of up to 
60 cm while the MSGC covers a radial region of up to 1.15 m. Again one sees a marked 
improvement over CDF in the region of interest. 

In addition, as can be seen from Fig. 23, the production cross section for a given mass is 
now increased by at least two orders of magnitude at the LHC because the center of mass 
energy is now 14 TeV. Such an increase in the cross section combined with the increase in 
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the radial distance covered by silicon detectors would, in principle, help in the search for 
long-lived quarks. One might wonder if a similar analysis as the one preformed by CDF 
could be carried over to the LHC experiments. Considering the fact that, with a much 
higher energy and consequently larger cross section, the number of events and background 
will be significantly higher as well. This would probably require a different search algorithm. 

Finally, concerning proposed but not yet approved colliders such as the Next Linear 
Collider (NLC) with y/s = 500 GeV, the long-lived quarks with the mass range discussed in 
Section IV, would be produced copiously and with little background. What kind of signal 
would one search for will depend on the kind of detectors involved. Whether or not the 
existence (or nonexistence thereof) of these long lived quarks will be established by CDF, 
ATLAS, or CMS by the time the NLC operates (if approved) remains an open question. 

B. Lepton Searches 

Earlier in this Report, indirect bounds on the masses of heavy leptons arising from 
violations of e — /i — r universality were discussed. The bounds were very sensitive to the 
mixing angle between the third and fourth generations. In this section, we discuss direct 
detection of heavy leptons. 

All current experimental bounds on heavy leptons come from experiments at electron- 
positron colliders. This is not surprising; the cross section for heavy lepton production at 
hadron colliders is small and backgrounds are large. Of course, once LEP200 shuts down 
in a couple of years, the only available colliders for searching for heavy leptons will be the 
Tevatron and the LHC. We will first examine the current bounds on heavy lepton masses, 
and then turn towards the future. 

It is generally believed that charged heavy leptons can be excluded up to the approximate 
kinematic limit of LEP. This is not necessarily the case however; the charged heavy leptons of 
many of the most interesting models have not been excluded for masses above approximately 
45 GeV. Below 45 GeV, heavy leptons (charged or neutral) would contribute to the decay 
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width of the Z, and such leptons can be excluded. 

The strongest bounds on heavy lepton masses reported by LEP have been reported by 
OPAL ||276|| and by L3 ||277|| . Both experiments assume that the heavy leptons decay via 
the charged current decay — the decay E — > tZ, which can be large in vectorlike models, is 
not considered (since it seldom dominates the charged current decay, their bounds are not 
affected). In the OPAL analysis, they exclude charged leptons which decay via E — > v{W 
with masses below 80.2 GeV, and those which decay via E — > NW with masses below 81.5 
GeV. Unfortunately, this latter decay assumes that the mass difference between the E and 
the N is greater than 8.4 GeV. As we have seen in earlier sections, vectorlike models have 
mass splittings on the order of a few hundred MeV, and even in chiral models, the splitting 
could also be small. Note, however, that if the mixing angle between the third and fourth 
generation is bigger than 10~ 6 , then the E — > v{W would occur near the vertex, and the 
OPAL bound would apply. The bound was obtained for LEP at \/s = 170 — 172 GeV, and 
can be improved somewhat for the later runs. 

In the L3 analysis, the mass splitting was assumed to be larger than in the OPAL case, 
greater than 10 GeV, and similar bounds were obtained. The L3 analysis also looked for 
long-lived charged leptons, which would exist if the mixing angles with lighter generations 
were small (typically less than 10~ 7 ) and the charged lepton is lighter than its neutral 
partner. Of course, such leptons must eventually decay, for cosmological reasons, but we 
discussed a variety of such scenarios earlier in this Report. L3 excludes such leptons up to 
a mass of 84.2 GeV. 

Both experiments also looked for heavy neutrinos, which decay at the vertex (mixing 
angle greater than 10 -6 or so) into a charged lepton and a W. For both experiments, 
the bounds for Dirac (Majorana) neutrinos are approximately 78 (66) GeV for decays into 
electrons or muons and 70 (58) for decays into taus. 

So, summarizing the current situation, the bounds on the charged heavy lepton are 
approximately at the kinematic limit of the collider if and only if this lepton is either stable 
(i.e. with a lifetime greater than tens of nanoseconds), has a large (8 GeV or greater) 
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splitting with its neutrino partner, or has a relatively large mixing angle (10~ 6 or greater) 
with lighter generations. Note that one of the most interesting models is the E 6 motivated 
model with a vectorlike doublet with very small mixing, and this lepton satisfies none of 
the above conditions. The mass bound for such a lepton is still only given by Z decays. 
(A search for a nearly degenerate lepton doublet was reported [ |278| ] many years ago by the 
Mark II detector, but only applied for leptons lighter than 10 GeV.) 

Could more analysis at LEP find such a charged heavy lepton? In vectorlike models, the 
principal decay of the E is into the N plus a very soft pion. It appears to be impossible 
to pick this pion out of the background from soft tracks from beam-beam interactions. 
Recently, Thomas and Wells ||181|| proposed a new signature — triggering on an associated 



hard radiated photon. This is similar to proposals for counting neutrino species through 
e + e~ — > 1/P7. At LEP, one would look for e + e~ — > L + L~^y. There are backgrounds from 
the above neutrino process, but they can be reduced by looking for a displaced vertex (the 
decay length is of the oorder of centimeters) and for the soft pions. Thomas and Wells plot 
the cross section as a function of the L mass and the minimum photon energy. With an 
integrated luminosity of 240 pb _1 at y/s = 183 GeV, and a minimum photon energy of 8 
GeV, they estimate that a doublet mass of up to 70 GeV could be detected. Presumably, 
this reach will be considerably higher for the more recent higher energy runs. They note 
that this signature is unusual for e + e~ machines since it is not only limited by the machine 
energy, but also by the luminosity, and that higher luminosity can significantly extend the 
reach. 

In the chiral case, there remains a "hole" to be filled. If the E is either very close in 
mass to or lighter than the N, it will primarily decay via mixing. If the mixing angle is 
greater than about 10~ 6 , then it decays near the vertex and can be detected at LEP up to 
the kinematic limit. If the mixing angle is smaller than about 10~ 7 , it is effectively stable 
and can be detected at LEP up to the kinematic limit. For intermediate angles, the decay 
length is of the order of tens of centimeters to a meter. Of course, some will still decay near 
the vertex, and some will decay well within or outside the detector, and so it is possible 
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that a complete analysis could close this hole. Doing so would be useful, since one of the 
plausible values for the mixing angle discussed in Section III is given by ^Jm UT /M N , which, 
using 0.01 eV for the v T mass, gives 3 x 10~ 7 for the mixing angle. 

In the future, similar analyses to the above will give similar bounds at lepton colliders 
near the kinematic limit of the colliders. We have noted, however, that the detection of 
the vectorlike doublet leptons remains problematic and can best be attacked looking for an 
associated hard photon. There also seems to be a window for decay lengths of the order of 
tens of centimeters which has yet to be closed. 

The decay mode E —>■ tZ is generally smaller than the charged current decay-although 
it is a much cleaner mode, backgrounds are not the problem for E — > uW, thus the latter 
would be detected first. This neutral current decay mode, however, provides a much cleaner 
signature for hadron colliders, which we now discuss. 

A study of searches for heavy charged leptons at hadron colliders was performed by 
Frampton et al. ||182|| They considered charged lepton production at the SSC and at the 
LHC (at 17 TeV). There are two main productin mechanisms for heavy leptons at hadron 
colliders. The first is gluon fusion, through a triangle graph, into a Higgs boson or a Z-boson. 
The second is quark fusion directly into a Z-boson (the effects of photon exchange are much 
smaller than those of the Z). The cross section for lepton production through quark fusion 
falls off very rapidly as the lepton mass increases, but the cross section through gluon fusion 
does not fall off as rapidly, since the matrix elements increases as the square of the lepton 
mass. 

First consider the chiral case. Here, gluon fusion dominates for lepton masses above about 
150 GeV, and the total cross section for masses between 100 and 800 GeV drops from 0.5 pb 
to 0.05 pb. This will lead to many thousands of events per year at the LHC. The signature 
would be a conventional heavy lepton signature. For those masses, and for chiral leptons, 
one can expect a reasonably large splitting between the N and the E, leading to standard 
single lepton and missing momenta signatures; even if there was an unexpected degeneracy 
(or if the N were heavier) mixing would lead to clear signatures (note, as discussed above, 
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the importance of closing the window for mixing angles near 10~ 6 ). 

What about the vectorlike case? Here, gluon fusion doesn't contribute, since the leptons 
don't couple to the Higgs and the vectorlike coupling to the Z gives no contribution due to 
Furry's theorem. Thus, the contributions are only through quark fusion, which fall off much 
faster. As the lepton mass increases from 100 GeV to 800 GeV, the cross section falls from 
1 pb to 0.001 pb. For a 400 GeV heavy lepton, this will give only 1000 events annually at 
the LHC. This makes detection more difficult, however one should recall that these leptons 
can decay via the neutral current: E — > rZ which has a branching fraction of at least a 
few percent (and in some models much larger). This would give a very clear signature with 
very low background. Even if the decay is suppressed by very small mixing angles, and thus 
the E passes through the detector, stable lepton searches should see it (it can be readily 
distinguished from a muon by time-of-flight using the velocity distributions as given in Ref. 

EH-)- 

VIII. CONCLUSIONS. 

There are still several reasons to believe that further quarks and leptons remain to be 
discovered. Although the fourth or further quark-lepton generations cannot be exactly sim- 
ilar and sequential to the first three generations, there are plenty of alternative possibilities 
which avoid the experimental embarassment to the invisible Z partial width of a fourth light 
neutrino. The additional quarks and leptons may be chiral as in the first three generastions 
or non-chiral and vector-like. 

The allowed masses are constrained by the precise electroweak data particularly at the Z 
pole where the data now agree with the minimal standard model at an astonishing 0.1% level. 
The S, T, U parameters then restrict what states may be added as discussed above in Section 
(III). Also the stability of the observed vacuum places constraints on additional fermions as 
does the (optional) requirement of grand unification of the three gauge couplings. Mixing 
angles for the new quarks and leptons are relatively unconstrained, except by unitarity, 
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without new experimental data. 

The lifetime and decay modes (see Section (IV)) of a heavy lepton depend critically on 
whether the N or E state is the more massive. A similar dependence occurs for heavy quarks 
which may have such small mixing with the known quarks that at least one new quark may 
have an exceptionally long lifetime. 

In Section (V) we have considered the fascinating possibility that the Higgs boson is 
not elementary but rather some bound state of additional fermions which transform under 
the standard gauge group. The heaviness of the top quark has suggested to some that it 
plays a special role in electroweak symmetry breaking, but even heavier fermions are more 
attractive candidates to participate in dynamical symmetry breaking. 

CP symmetry violation has two disparate but likely related aspects in the Standard 
Model: the strong CP problem and the weak CP violation in kaon decay. Strong CP can be 
addressed by addition of extra quarks as explained in our Section (VI). Weak CP violation by 
the KM mechanism requires at least three generations, and acquires even more CP violating 
phases in the presence of additional quarks. We have illustrated this with the Aspon model 
which invokes spontaneous CP violation to relate solution of the strong CP problem by 
extra vector-like quarks to the violation of CP symmetry in kaon decay. The new vector-like 
quarks may have long lifetime as mentioned in Section (IV). 

Experiment is the final arbitor of everything we have reviewed. Long-lived quarks are 
being sought at collider facilities. To some extent, detectors have not been designed for such 
a possibility and this review may encourage further thought in detector design. Similarly 
heavy leptons are being, and will be, investigated at existing and future colliders. 

Discovery of a further quark or lepton would be revolutionary and propel high-energy 
physics in a new and exciting direction. 

This work was supported in part by the US Department of Energy under Grants No. 
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